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Abstract 

We provide a detailed analysis of Friedmann-Robcrtson- Walker universes in a wide range 
of scalar-tensor theories of gravity. We apply solution-generating methods to three 
parametrised classes of scalar-tensor theory which lead naturally to general relativity 
in the weak- field limit. We restrict the parameters which specify these theories by the 
requirements imposed by the weak-field tests of gravitation theories in the solar system 
and by the requirement that viable cosmological solutions be obtained. We construct a 
range of exact solutions for open, closed, and flat isotropic universes containing matter 
with equation of state P < and in vacuum. We study the range of early and late-time 
behaviours displayed, examine when there is a 'bounce' at early times, and expansion 
maxima in closed models. 



PACS number(s) 98.80.Cq, 98.80.Hw, 04.50.-hh, 04.60.-m, 12.10.-g 



I Introduction 



Cosmological models arising from theories of gravity in which the Newtonian gravi- 
tational 'constant', G, varies with time have a long history. They were first studied in 
detail in response to Dirac's claims that a coincidence between the values of 'large num- 
bers' arising in dimensionless combinations of physical and cosmological constants could 
emerge naturally if one of the constants involved possessed a time- variation that was signif- 
icant over cosmological timescales |]l[|. Dirac ascribed that time variation to G and simply 
wrote the time-variation into the Newtonian expressions which held for constant G. Subse- 
quently, mathematically well-posed gravitation theories were developed in which Einstein's 
theory of general relativity (GR) was generalised to include a varying-G by deriving it 
from a scalar field satisfying a conservation equation. These scalar-tensor gravity theories, 
first formulated by Jordan Q, were most fully exploited by Brans and Dicke in 1961 
Motivated by claims that the observations of light-bending by the sun were in significant 
disagreement with the predictions of GR, Brans and Dicke explored the possibility that 
the simplest scalar-tensor theory could provide predictions of the weak-field solar-system 
tests in agreement with light-bending and other geological and paleontological observations 
IQ]. Cosmological models could be found in Brans-Dicke (BD) theory, but astronomical 
observations were unable to impose stronger limits upon them than had been found from 
solar system experiments. Subsequently, the doubts regarding the compatibility between 
observations of solar light-bending and the predictions of GR were removed by a fuller un- 
derstanding of the uncertainties surrounding measurements of the solar diameter at times of 
high solar surface activity |^] . This removed the one observation that called for the replace- 
ment of GR by a scalar-tensor theory displaying varying-G in the weak- field limit. Since 
that time there have been two observations which have led to renewed interest in gravity 
theories with non-Newtonian variation in G: one was the claim that a 'fifth force' of Nature 
existed on laboratory length scales and influenced Eotvos experiments as if there existed 
a deviation from the inverse-square form of the law of gravitation in the non-relativistic 
Newtonian regime the other has been the occasional claim that the flatness of the ro- 
tation curves displayed by spiral galaxies might be a signal of non-Newtonian gravitational 
attraction |^. The evidence for the 'fifth force' variations has not been supported by other 
experiments and flat rotation curves appear as natural outcomes of protogalaxy formation 
in conventional gravitation theories. In both cases, the deviations from conventional gravi- 
tation theory, with constant G, would introduce a new characteristic length scale into the 
law of gravity with no fundamental basis other than to explain a particular set of obser- 
vations. Some attempts were made by Gibbons and Whiting to see how simple scaling 
arguments might relate preferred high-energy physics scales to those observed in fifth force 
experiments. More recently, the study of scalar-tensor gravity theories has been rejuvenated 
by theoretical developments in the study of the evolutionary possibilities open to the early 
universe. 

In a metric scalar-tensor theory of gravity the gravitational coupling is derived from 
some scalar field, (p, so G = G{(f)). Historically, most interest has been focused upon the first 
and simplest theory of this type, presented by Brans and Dicke [^, in which the coupling 
function a;((/>) is a constant. In general, the consideration of scalar-tensor theories with 
non-constant a;(0), § greatly enlarges the range of possible G- variations and weakens 
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the impact of observational limits accordingly. In the weak-field limit Nordtvedt found an 
expression for the observed value of the gravitation 'constant' in these theories, to leading 



order, as |10|. 



G{t) = cp- 

so that 



2uj( 



3 + 2lj( 



G ■ ('i + 2uj\ ( 2w'(0) 
G + 



G V4 + 2u;; V (3 + 2(j)2 

We see that in Brans-Dicke theory, with constant uj the variation of G{t) is just inversely 
proportional to Moreover, for one particular choice of w(^), 



2(<A-1)' 

it is possible to have G = to first order in the weak-field limit . 

Our observational limits are a mixture of cosmological limits, studies of astrophysical 
objects, and weak-field tests of gravitation in the solar system. In Brans-Dicke theory the 
scope for significant deviations from constant G is very small because of the constancy of 
Lo. However, if u varies then it can increase with cosmic time in such a way that a; — > oo 
and Lo'oj~'^ ^ as t — > oo so that weak-field observations at the present time accord well 
with the predictions of GR even though the theory may possess significant deviations from 
GR predictions at very early cosmological times [IT^j . 

The principal observational upper limits on G jG come from many different observations. 
Passive radar data on the distances to Mercury and Venus 14 1 give the limit GjG < 



4 X 10"^°yr-^ This was improved using the Mariner 9 Mars orbiter by Anderson et al. 
1 15 1 to G/G < 1.5 X 10^^'^yr^^. Anderson et al. |lH] then used Mariner 10 data and radar 



ranging to Mercury and Venus to obtain G/G < 0.0 ± 2.0 x 10 ^"^yr ^. Viking landers. 
Mars orbiters and transponders gave upper limits of G/G < 3 x lO^^^yr^^, ll^], G/G < 



2 ± 4 X IQ-^^yr-^, |l|], and G/G < -2 ± 10 x lO^^^^^-i^ ^ xhese are made uncertain 



by our incomplete knowledge of the asteroid distribution. Studies of the Binary pulsar PSR 



1913+16 by Damour, Gibbons and Taylor |0[ give a limit of G/G = -(1.10 ± 1.07) x 
10~^^?/r~^ with uncertainties due to the pulsar's proper motion. Gravitational lensing 
promises to yield new tests, but depends upon other uncertain cosmological parameters 



|21|. Krauss and White argue that a variation in G of AG/G < 20 would infiuence lenses 
at z ~ 1.5 and these studies might eventually achieve constraints of order G/G < 10"^^ 
in Brans-Dicke theories. White dwarf cooling will be affected by variations in G, [^], and 
recently, Garcia-Berro et al. gave limits, depending upon the chemical composition 
of the white dwarf of G/G < —(1 it 1) x lO^^^yr^^ for chemically stratified models, and 
G/G < — (3]'^3) X 10~^^yr~^ for models with constant C/0 composition. However, these 
studies have theoretical uncertainties introduced by the quasi-Newtonian modelling of the 
variation of G, [24|. This is also true of attempts to constrain G/G by measurements of 
the time- variation of neutron star masses [25|. Cosmological nucleosynthesis gives limits 



that are of order G/G < O.OIH in Brans-Dicke theory. Limits have also been obtained in 
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theories with 2a; + 3 cx 



which display a slow logarithmic decrease in G{t), 

1 fjrt 



-10 



yr 



Other 

lunar occultations 
lunar laser-ranging studies of the Moon's 
< 0.3 X 10-^°2/r-\ and Miiller et al. [0] derive 
G/G < (0.1 lb 10) X 10~^^yr~^. There are proposals to measure via radar the position of a 
transponder in orbit around Mercury, [^]. The ranging sensitivity is claimed to be able to 
reach just a few centimetres which would translate into limits of order G/G < 3 x lO^^^yr^^. 
Most of these limits (especially cosmological ones from nucleosynthesis) are uncertain if 



limits include those from solar evolution, |28], of G/G < 10 
and eclipses, [||], of G/G < 0.4 x 10" 
orbit around the Earth, pi, oi G/G 



varying G is coupled to the variation of other constants |33]. 

The coupling function of scalar tensor theories, a;(0), is related to the PPN parameters 
by, lie 
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Observational limits on the weak field PPN parameter 7 are 7 ~ lit 0.002 from radio timing 
delays [p4| ; 7 ~ 1.0002ib0.002 from light deflection using VLBI observations of quasars ||3^ ; 
7 ~ lib0.02 from lunar laser ranging, |31] . Future experiments, GPB, POINTS and Mercury 
Relativity Satellite, hope to reach sensitivities of I7 — 1| ~ 3 x 10~^, ||3^. For combinations 
of two PPN parameters determining uj, limits of |4/? — 7 — 3| < 5 x lO^^(lcr) have been 

I7- 1| < 0.002 and |4/3 



found [19|. If we take the observational limits as 
then we have uj > 498 and 



7-3| < 0.001 



ijj 



(3 + 2a;)2(4 + 2L.;) 



< 0.001 



where u)' is evaluated at the asymptotic (f. 
a rather weak limit of \uj'\ < 10~^ x 999^ 



value 
-0(1) 



where 00 



bo). Hence, we have only 



Planetary data |37| also provides a limit on the spatial gradient of G over solar system 
scales of VG/G < 3 x 10~^°AC/~^. A limit of 6G/G < 10~^^ on possible spatial anisotropy 
of G in the solar system has been derived by studying the alignment of the Sun's rotation 



axis with the direction of the solar system's angular momentum vector |38| , and of 5G/G < 
2 X 10^^^ from satelhte and LAGEOS laser ranging data ||39t| . 

A brief period of interest in extended Kaluza-Klein cosmological theories ||4 



which 

culminated in their replacement by superstring theories [41|, revealed how time- variations 
in any extra (> 3) dimensions of space would manifest itself through the time- variation of the 
'constants' defined in the three dimensions we observe. Over the same period, cosmologists 
showed growing interest in the behaviour of all scalar fields during the early stages of the 
universe. The time variation of a field energy source in the early universe at a rate slower 
than the universe is expanding is a general possibility only for scalar fields and has become 
known as 'slow-rolling' of the field. Typically, it produces an acceleration of the expansion 
scale factor of the universe with time: a phenomenon known as 'inflation' [42|. This has 
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led to the investigation of general-relativistic cosmological models containing a wide range 
of self-interacting scalar- field sources [43|, the classification of different varieties of inflation 
that can result from their slow-rolling evolution, and the extraction of detailed predictions 
concerning the fluctuations imprinted in the cosmic microwave background radiation by 
the spectra of gravitational waves and density perturbations that emerge from a period of 
primordial inflation. However, it has also been recognised that natural scalar fields might 
be provided by the scalar component of a scalar-tensor theory of gravity |44]. Such theories 
possess close conformal relationships with GR plus explicit scalar fields, and although the 
scalar field determining the strength of the gravitational coupling does not easily drive 
inflation it does have a significant effect upon the pace of inflation that arises when self- 
interacting scalar flelds are included in the universal energy-momentum tensor. Following 
the long initial study of BD cosmological models. Barrow |12] showed how to generate 
cosmological solutions of vacuum or radiation-dominated isotropic cosmological models in 
any scalar-tensor gravity theory. The method used works for any isotropic cosmological 
model with an energy momentum tensor possessing a vanishing trace. Barrow and Mimoso 
1 27 1 then devised a more complicated procedure which allows solutions to be generated for 
isotropic, zero-curvature universes with the pressure p and density p related by a perfect- 
fluid equation of state p = (7 — l)p, with the constant 7 lying in the range < 7 < 4/3. 
In particular, they found the first dust solutions for scalar-tensor theories more general 
than BD together with a wide range of new inflationary solutions. Recently, the techniques 
introduced by Barrow and Mimoso have been used to study the asymptotic behaviour 
of isotropic and anisotropic cosmological models in scalar-tensor theories 46 1. The 
anisotropic case has also been studied by Mimoso and Wands |47]. In this paper we are 
going to extend these studies to arrive at a more general and systematic understanding 
of the behaviour of isotropic cosmological models containing matter with equation of state 
P ^ /o/3. We will be interested in studying scalar-tensor gravity theories which can approach 
general relativity in the weak-fleld limit at late cosmic epochs. By studying the classes of 
gravity theory which allow this approach to occur we shall use our solution-generating 
techniques to build up a detailed picture of the behaviour of scalar-tensor cosmological 
models. We shall be particularly interested in models containing radiation {p = p/3), dust 
(p = 0), and inflationary stresses {p = —p). 

The plan of the paper is as follows. In section |ll| we introduce the Lagrangians and fleld 
equations which define scalar-tensor gravity theories in terms of their coupling function 
ijj{4>) and dictate their evolution. We shall specialise our study to the cases of Friedmann 
universes containing perfect fluid sources. In section III we describe the two techniques for 
finding complete and asymptotic solutions of these equations for arbitrary choices of uj{4'). 
The 'direct' method works for universes of all curvatures but is restricted to the vacuum and 
radiation-dominated cases. The 'indirect' method works only for flat universes but for any 
equation of state. We shall be especially interested in the vacuum, dust (p = 0), radiation 
and inflationary {p = —p) cases. These techniques allow us to draw important conclusions 
about the small and late-time evolution of cosmological models in scalar-tensor theories. 
In section IV we introduce three broad classes of gravity theory, defined by the form of 
Lo{(j))^ which allow distinct forms of cosmological evolution. The free parameters can be 
restricted to allow the theories to reproduce the successful weak-field predictions of general 
relativity and to give cosmological solutions. The cosmological consequences of these three 
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classes of theory are explored systematically in sections 0, [V^ , and VD . In each case we are 
interested in determining the early and late time behaviours, and finding exact solutions 
which describe the dust, radiation, and inflationary phases of expansion. The results are 
discussed in section IVID . 



II Scalar-tensor cosmologies 

A Field equations 

We shall consider scalar-tensor theories of gravity defined by the action, 

(1) 



Sg = I d^Xy 



AT? _L abo 

-(/)7C H —g da 

<p 

where TZ is the curvature scalar arising from the spacetime metric ^(ab, 9 is the determinant 
of (7ab, (j) IS a scalar field, and u;{(p) is a function determining the strength of the coupling 
between the scalar field and gravity. We are working in units such that Newton's constant. 
Gat, is equal to unity. The field (f> is the analogue of G~j^^ in GR except that here, in contrast 
to Einstein's theory, (pis a dynamical quantity. Scalar-tensor gravitational theories therefore 
permit histories in which the value of the gravitational "constant" varies. The simplest such 
case is that explored by Brans and Dicke, where w is a constant. The action in Eq. (|l|) 
offers more general Lv{(j)) theories as natural extensions to BD gravity. It is these theories 
that will be of primary interest in this paper. Demanding that the first-order variations of 
Eq. (H) with respect to (/) and g^h vanish, we derive the field equations: 



^ab - ^gahT^ = -Svr^ - (^0a(Ab " ^^ab^c'?^'''^ " ^ (0a;b " 9ah^(l>) , (2) 

and the conservation law, 

r% = o, (4) 

where T^ab is the Ricci curvature tensor. Tab is the energy-momentum tensor specifying the 
properties of the matter occupying the universe. Primes indicate derivatives with respect 
to (p. The first of these three equations is the scalar-tensor analogue of Einstein's equations, 
the second is the wave equation for 0, and the final expression is the energy-momentum 
conservation law for the matter, which ensures that each theory is consistent with the 
equivalence principle. 



B Friedmann universes 



We shall examine solutions to these equations that describe homogeneous, isotropic 
Friedmann-Robertson- Walker (FRW) cosmological models, with time-varying-G. The FRW 
metric line element in spherical polar co-ordinates (t, r, 6, ip) is given by: 



dr'' 
1 — kr^ 



+ r' 



(de"^ + sir^edi}'- 



(5) 
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where the curvature parameter is k = —1,: 0,: +1 for open, flat or dosed cosmologies re- 
spectively, and the scale-factor, a{t), characterises the expansion history of the universe. We 
shall assume the universe contains a simple perfect-fluid which may be accurately described 
by a perfect-fluid equation of state. 



p = (7 — ; < 7 < 2, : 7constant ; 



with these prescriptions, the equations of motion become. 



k 



Svr p 
3" 



(6) 



(7) 



+ 



3H + 



UJ( 



2oj{(l))+3 



2uj{(j)) + 3 



(4-37), 



H 



<j) Snp (37 - 2)uj{(f)) + 3 1 uj{(p) 



3(j) 2uj{(j)) + 3 



+ 



2 2uj{(p) + 3 < 



(8) 
(9) 



p + SjHp = , 



(10) 



where H = a/a is the Hubble expansion parameter and overdots denote derivatives with 
respect to comoving proper time, t. 

The structure of the solutions to these equations is sensitive to the form of the coupling 
function w (</>), which deflnes the scalar-tensor theory of gravity. As Nordvedt flrst showed 
1 10 1, it is possible to place bounds on the parameter-space of many models, prior to extract- 
ing a full solution from the field equations, simply by inspecting the explicit form of (^{(p). 
Nordvedt 's constraints demand that the theories tend to GR in the weak-field limit, so that 
they concur with the observational limits on light-bending and perihelion precession. This 
requirement manifests itself explicitly in the conditions a; — > 00 and uj'uj~^ — > as t — > 00. 
Typically, scalar tensor theories add a term proportional to lo'uj~^ to the weak-field predic- 
tions of general relativity. Whilst the first condition (cj ^ 00) is well known, the second 
{u!'uj~'^ 0) is not. As we shall see in Section IV there are many theories which satisfy 
the first condition but not the second: for example those defined by uj{(j)) ~ (1 — i;^>/0o)~" > 
with < a < 1/2, as (/> — > (/)o from below. 



Ill Methods of solution 



Solutions to Eqs. (0) - ( [lO|) for Brans-Dicke FRW models have existed for some years [48, 
49 1 . The Brans-Dicke case has also been studied qualitatively by Kolitch and Eardley, |50|. 



Recently, Barrow ||T^, 52 1 and later Barrow and Mimoso have extended these treatments 
to general (jo{4>) theories, providing a method for obtaining zero-curvature cosmological 
solutions when 7 < 4/3, and solutions for vacuum and radiation-dominated universes of 
any curvature. 

We now recapitulate the two methods of solution introduced in refs. [^] and |27| to 
solve the cosmological field equations for theories specified by any uj{4>)- 
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A Vacuum and radiation models 



1 Exact solutions 

The general solutions to Eqs. - ( [lO| ) contain four arbitrary integration constants, one 
more than their GR counterparts, the extra degree of freedom being attached to the value 
of (p. When the energy-momentum tensor is trace-free there exists a conformal equivalence 
between the theory and GR, the right-hand side of Eq. (^) vanishes and = is always a 
particular solution, corresponding to a special choice of the additional constant possessed 
by the model over GR. Consequently, the exact general solution of Einstein's equations 
when Tab is trace-free is also a particular solution to Eqs. (0) - (|To| ) with and hence a;((/>), 
constant. 

It will seldom be the case that the particular solution obtained in this way will form the 
general solution for that particular choice of uj{(I)) ■ Usually, however, it will be the late or 
early time attractor of the general solution. For example, in the case of Brans-Dicke theory 
the special GR solution is the late-time attractor for flat and open universes but not the 
early-time attractor. However, one of these authors |12| developed a method for integrating 



the field equations for models containing trace-free matter. The procedure is as follows. 
Eq. ( p!o| ) integrates immediately to yield 

87r/9 = STa"^^ , (11) 

where F > is a constant of integration; F = describes vacuum models. Making the 
choice 7 = 4/3, corresponding to blackbody radiation, and introducing the conformal time 
co-ordinate, rj, defined by 

adri = dt , (12) 

Eq. (P) becomes 

where subscript r] denotes a derivative with respect to conformal time. This integrates 
exactly to give, 

(P^a"^ = 3'^/'^A{2io{(l)) + 3)~^/2 ; A const. (14) 



We now employ the variable used by Lorenz-Petzold to study Brans-Dicke models [^], 

y = (pa\ (15) 
to re- write the scalar-tensor version of the Friedmann equation, Eq. (0), as 

{y^f = -Aky' + 4Fy + ^ {(t>.,f a\2iv{<P) + 3) . (16) 



Dividing Eq. (|T^) by Eq. (|T^), and using Eq. (14), we obtain the coupled pair of differential 
equations 

% = 3V2Ay-l(2^(<A)+3)-l/^ (17) 
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We may now obtain the general solution for a particular choice of uj{(j)), given k. Integrating 
Eq. ( |l8|) yields y{r]) which, in conjunction with uj{(I)), implies 4>{r]) from Eq. ( p^ and, without 



further integration, a{ri) from Eq. (15 ). If Eq. (|T2|) is both integratable and invertible we 
may compute and a{t), so completing the solution. 

The vacuum models are obtained by setting F = and in this case Eq. ( p^ ) has three 
possible solutions, according to the value of k: 



yin) 



( A{r] + r]o), k = 0, 

iAsinh[2(r/ + 7?o)], k = -1 , 
^ iAsin[2(r/ + r/o)], A; = +1 , 



(19) 



where tjq is an arbitrary constant fixing the origin of ij-time. Combining these results with 
Eq. ([TtI) yields the set of integral relations: 



(2cj((/.) + 3)V2 



\/31n|7? + 7/o| , k = 0, 

\/3 In |tanh(r/ + ryo)! j k = —1 , 
\/3 In |tan(r/ + r/o)| ; k = +1 . 



(20) 



Specification of uj{(f)) allows the full solutions to be completed. When k = —1 the negativity 
of the right-hand side, arising because < | tanh(r/ + ?7o)| < 1 V r/, places strong constraints 
on the allowed form of the integral on the left. 

Similarly, one may perform this treatment on the radiation models, ie those cases for 
which r > 0. Again, the results are classified by the value of k: 



( r(r? + Tjo)^ - Ay AT . 



k = 0. 



ir + i(A2 - r2)i/2 sinh[2(7? + r?o)] , k = -1 



(21) 



^ ir+l(r2 + ^2)1/2 gi^[2(r? + 77o)]. 



Integrating Eq. (jl^) with the solutions above leads to 

\/31n 

(2u;(<^)+ 3)1/2 



2r??+2rr?o-A 
2r»7+2rr?o+A 



rtanh(??+>?o) + (A^-r^)l/^-yl 
rtanh(»y+»7o)+{A2-r2)l/2+yl 

rtan(r;+>;o)+(r2+A2)V2-A 
rtan{r]+rio)+{r2+A'^y/^+A 



+ 1. 

k = 0, 
k = -l, 
k = +l, 



(22) 



These expressions can be simplified by choosing the arbitrary integration constant r]o such 
that 2r??o = A. 

The domain on which each right-hand side exists strongly constrains the integral on the 
left. For instance, if we require (j) £ (0,<;/'o) then the corresponding range of the function 
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of (p resulting from the integral on the left must be compatible with the allowed range on 
the right. In general, (j) must tend to its general relativistic form, ie a constant, at large 
rj. We shall find in Sections ^ - VIl that this behaviour is not generic and often requires 



the integration constant associated with the left-hand integral to assume a particular value. 
This integration constant can be interpreted as an initial boundary condition on or w at, 
say, T] = 0. Further restrictions can be found by studying the evolution at early times. 
In Sections ^ - |VII we shall exploit these relations to derive general u;((/>) solutions for 



all values of the curvature parameter. 

2 Approximation techniques 

Many of the models we shall present are insoluble in terms of t. The reason for this is 
the non-invertibility of t{r]), arising from integrating a(r/). In these cases we shall invert 
t{ri) approximately at early and late times to obtain series solutions for the behaviour of 
(f){t) and a{t), indicating their limiting forms and their approach to these forms to leading 
order. We shall use the inversion technique of Olver [p3||. If we have an expression 



yix) = f{x)g{x) , (23) 

and we require an approximation to x{y), valid in a region where g dominates over /, we 
neglect / and write 

x{y)c^g-\y). (24) 

The next-order approximation is easily obtained by substituting this result into /, yielding 
/ ~ f{g~^{y))- Using this in Eq. (23) we have 



This procedure can be iterated indefinitely but we shall use the 2'^'^ iteration form as given 



in Eq. (H). 

To analyse the asymptotic behaviour we need to establish the value of r/ as t ^ oo and 
our primary interest is in models which display GR, ie (p ^ (j)o, in this limit. Radiation- 
dominated, A; = 0, universes in GR evolve like a oc t^^"^ at late times and hence t on rf and 
r] ^ oo as t ^ CO. There does not exist a spatially flat vacuum model in GR, however, 
the negatively curved general relativistic models are asymptotically vacuum-dominated, (as 
the matter density redshifts to zero) tending to the Milne solution, a (x t, as t — > cxd, and 
r/ oc Int which diverges with t. We therefore study the asymptotic behaviour of vacuum 
and radiation models in the limit r/ ^ oo when k < 0. We will not consider the late-time 
limit of the k > models. Typically, they recollapse to a final singularity if p > and 
p + 3p > 0, although a bounce occurs for many choices of uj{(j)). The behaviour in their 
recollapse phase is similar to the time-reverse of the early expansion of k = models. 

We now define the early-time limit. We examine first models with < (p < (po, in this 
range \/2ui + 3 does not change sign, as guaranteed by the choices of the coupling function. 



u}{<p), detailed in Section IV and Eq. (|14| ) ensures that (p is monotonically increasing. We 
thus extrapolate the evolution back to (p = and treat this as the early-time limit in 
vacuum and radiation models. In more general perfect-fluid models we find that we are 
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easily able to examine the behaviour when (j) > (j)Q hy exploiting the conformal invariance 
of the theory (this procedure is explained in detail in the Appendix) and we so examine 
early-time behaviour in the neighbourhood of the last zero, or non-zero minimum, of a. 



3 Maxima and minima 

It is easy to study the structure of non-singular vacuum and radiation-dominated cos- 



mological models within the framework presented in this section. Differentiating Eq. (15) 
and substituting from Eqs. (14) and (111) leads to 



{a\ = ^ - . (26) 

The condition for the scale-factor to contain a stationary point, = , is equivalent to 
{a^)n = when a ^ 0, which leads to the simple relation 

when (p is finite. This relation will prove useful for locating minima when we lack an exact 
solution for a(r/). 



Eq. (27) is a simple test for the existence of stationary points in the evolution of a. 



Stationary points of the type max'iSmm '^i^^ occur when a^,;<0, respectively but, because 
a > 0, it is sufficient to replace ajjrj with (a^).^^ in this condition. Differentiating Eq. ( p^ 
twice and evaluating at a stationary point, ie where (a^)^ vanishes, we obtain 



A condition for the stationary point to be a maximum then 

Uriri > 4'ririO, 

Differentiating Eq. (|l^) and discarding {a^)r) we have 



maximum 

(29) 



a2 (2a; + 3)" 

Substituting this into Eq. (p9|) , and remembering > yields 



(30) 



>_J^4V(0L (31) 
^''^ a2(2a; + 3)2' ^ ^ 



as the condition for the stationary point to be a max\mum- When < (p < (pQ we have 
uj'{(p) > for the choices in Section |^ and the right-hand side of Eq. (31) is negative 
definite. Thus whenever we can prove y^^ > we may exclude the possibility of a{rj) 
containing maxima, and then by continuity we can limit the number of minima to one. 
We test this condition for the forms of y{r]) given in Eqs. (H) and (|l]), the results are 
summarised in Table |^. From the table it can be seen that all spatially flat and negatively 
curved models may only contain a single minimum. The k = +1 models are not bounded 
in this way and may contain an undetermined number of minima and maxima. 
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Matter source 


Vacuum 


Radiation 


k 


-1 





+1 


-1 





+1 




> 


= 


< 


> 2r > 


= 2r > 


< 2r 



Table 1: UriTj for vacuum and radiation models of all curvatures. 

B General perfect-fluid cosmologies 

When T is non- vanishing the situation is substantially more complicated. In this in- 
stance, (/> = is no longer a particular solution of the field equations, forcing us to resort to 



more elaborate methods to obtain solutions. Barrow and Mimoso |27| have done this, for 



the = models, by generalising the method of Gurevich et al. [48| for BD models to the 
case of varying to. We now outline this procedure. 

Introducing the new time co-ordinate and the two new variables x and v such that 



X 



3(2-7), 



and confining attention to the k = Q models, Eqs. (0), (|8|) and (^) transform to 



2 ^ 2 
-x + v 



/2uj + 3 
V 3 



;2 + 4rA«3{2-7) 



v^ = T{4- 37) , 



and 



x^ = 3T [(2-7V + l] + ^( 



3{2lo + 3) 



X + V j LO^ , 



(32) 

(33) 
(34) 

(35) 
(36) 

(37) 



where subscript-.^ represents a derivative with respect to ^-time. Eqs. ( pq) and (37) integrate 
easily to yield 



r(4 - 37) (e - 6) 

3 



-V + V2UJ + 3 \^C + T{2--f) J V2uj + 3d^ 

C is an integration constant and ^1 fixes the origin of ^-time. Noting the relation 

3 



ad) 0^ ' " a 



nm) ^: — 



and differentiating y, with respect to ^, yields 

-MO 



— I + 



37-4 1 
— + 



(38) 
(39) 

(40) 

(41) 
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where a new function /{(,), is defined by 
r« 3(2-7) 



2u;((^) + 3 



/g, 2r(4-37) 
Solving Eq. we have the solution 



C + r(2-7) 



2a;(<?i) + 3c?C 



(42) 



6 5(0 



with g(^) simply related to f{^) by 



(43) 



(44) 



where D is a constant of integration. Eq. ( |40| ) immediately reveals a simple formula for the 
scale-factor: 



5(0^/(0 + ^^(e-6)' + ^: 



3 / 5 \ 



o = an - 



i ao constant . 



Finally, the scalar-tensor coupling function uj((j)) is given as a function of / by 

4-37 ur 



3(2-7)2 



/ + 



4-37 
3(2-7)^ 



/o 



(45) 



(46) 



where /o is another arbitrary constant. 

In the calculations to follow, we shall exploit the arbitrariness of and set it to zero. In 
addition, we note some relations between the other constants in the solution, arising from 
the scalar-tensor Friedmann constraint, Eq. (^5|) 



D 



rh ^3(2-7) 
00 «o 



4-37 
3(2-7)2 
r (4 - 37) , 



/o ) 



(47) 
(48) 



and from the requirement that the BD theory be recovered when uj = ujq is a. constant we 
obtain the further condition, 



,2_ /2r(4-37) 
V 3(2-7) 



/o 



(49) 



We shall require 0o, ie today, to be positive and from Eq. ( ^7| ) we see that this requires 
7 < 4/3. We shall also confine our attention to theories with 2a; + 3 > 0. Using Eq. ( p3[ ) 
we can re- write Eq. (^6|), as a function of (p and its derivatives, as 



2w(e) + 3 



4-37 



+ 



6-37 



3(2 - 7)2 (e A _ 30^^^ 



(50) 
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Asymptotically, 2u} + 3 ^ (55)^/5 as ^ ^ 00 and the above relation becomes 

2 

(51) 




The choice of </>(^) is thus equally fundamental as that of g{(,), and amounts to specifying 
G(^). We shall use the former. Once has been specified, we may infer g{^) from 



Eq. (|3|), a(^) from Eq. (|5|) and from Eq. (|^). If Eq. (|3|) can be integrated and 
inverted to yield ^(t) these variables may be expressed in terms of cosmic time, t. 

An important benchmark is provided by the behaviour of the BD theory, where uj{(j)) = 
ujQ = constant. In this case, the generating function, /(^), is given by a quadratic in ^: 



(52) 



Hence, in general (C / 7^ F) when 7 ^ 4/3, 2, we see that /bd oc as ^ ^ cxd and 
/bd oc ^ as ^ ^0, where dt oc a'^^"'~^^dS,. If we choose C = then /bd oc F^^ as ^ ^ 0. 
The choice C = restricts the solution to the special 'matter-dominated' solutions (termed 
'Machian' by Dicke js], ^, see also Weinberg |51]) which were first found for all perfect-fluids 



by Nariai [^]. If C 7^ then the early-time behaviour is dominated by the dynamics of the 
0- field; such solutions are termed '(/>-dominated' (or 'non-Machian' by Dicke). 

Therefore if we choose a generating function g{^) that grows slower than as ^ — > 00 it 
will produce a theory that approaches BD at late times ((/> ^constant, u;((/>) constant), 
whilst if g{(^) decreases slower than ^ as ^ — > then the theory will approach the behaviour 
of (/)-dominated BD theory at early-times. This means that we will find new (non-BD) 
late-time behaviours by studying generating functions which increase faster than g{^) = 
as ^ ^ 00 and new (non-BD) early-time behaviour by picking generating functions which 
decrease slower than g[^) = ^as^^Oor^^ ^min (if there is no zero of ^ at the minimum 
of a{t)). 



IV The coupling function 



We are interested in ascertaining the general behaviours displayed by cosmological mod- 
els in the range of scalar-tensor gravity theories that approach GR in the weak-field, late- 
time limit. This requires 2ti;-|-3— >ooast— >oo and also w'o;"^ — > if the solar system tests 
are to accord with observation. The specific form of the leading-order corrections to the 
general relativistic predictions of light-bending, perihelion precession, and radar echo delay 
are all almost equal to u!'(jJ~^ in the large uj limit. However, the rate at which tends 
to infinity will determine the form of the cosmological models. In an earlier paper we 
explored the behaviour of simple power-law forms for uj{(j)) which, although growing with 
time, only attain the GR limit when (p = 00, although at any finite time lo can be made as 
large as we wish by the choice of the constants defining io{(p). Here, we turn our attention 
to a potentially more interesting situation in which oj ^ 00 as (j) ^ (po where 4>q may be 
taken as the present value of (j){t), which determines the observed value of the Newtonian 
gravitation constant, G = (pQ^ . 
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Theory 1: 2cj(0) + 3 = 2Si|l - 



a) 



b) 



Figure 1: Theory 1 with a.) Bi = 5 and a = 0.5, 1, 1.5; and b) a = 2 and Bi = 5, 10, 15. 



We shah study the three general classes of theory, some examples of which are displayed 
in Figs. |l|-|: 



Theory 1. 2uj{(j)) + 3 = 2Bi\l 



'o| °'; a > 0, Bi > constants. 



Theory 2. 2a;(</>) + 3 = B2 |ln(</>/0o)r^^ 6 > 0, B2 > constants. 
Theory 3. 2Lj{(j)) + 3 = B3 1 - (0/(/)o)^ ~\ /3 > 0, > constants. 



Theory 1 has been studied previously by Serna and Alimi [^], Comer et al. |46], Barrow 
52|, and Garcia-Bellido and Quiros [^]. Serna and Alimi ^5| paid particular attention to 
the radiation eras of these models, extending their treatment to allow 2a; + 3 < 0, leading to 
(/)-dominated initial conditions. Their analysis dealt mainly with the early-time behaviour 
of the models, we shall examine in detail the late-time approach to GR. Theories 1-3 will all 
permit u; — > 00 if ^ c^q at late times and span a wide range of different rates of approach 
to GR in the weak-field limit. They can all be reduced to Barker's constant-G theory |ll[| 
as (p ^ (pQ for special parameter choices. Theories 1 and 3 approach Brans-Dicke theory 
1^] as (/> ^ and the power law form studied in |12| as ^ 00. General functional forms 
for 2ij((p) + 3 can be expanded in a series of functions of this form and their asymptotic 
behaviours at small and large times will be dominated by one term of the above type. In 
fact, further restrictions can be placed upon the allowed theories within these three classes 
by the weak-field limit requirements, as follows: 

Theory 1: We see that w — > 00 as ^ 0o is guaranteed if a > and that 



to 



^ oc 1 - — 



2a-l 



(53) 



Hence lo'uj ^ — > and the weak field limit will be compatible with observation as 
so long as a > 1/2. 
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Theory 2: Here, lv —>■ oo as (p —>■ (po for 5 > 0, but 



and this tends to zero as ^ so long as (5 > 1/4. 
Theory 3: Here, lo ^ oo as 4> ^ 4>o, but 



(54) 



UJ 



oc 



1 



<t>oJ 



13-1 



(55) 



and this tends to zero for all /? as — > (/>o- 

These constraints on the parameter ranges of theories 1-3 are independent of the form 
of the cosmological solutions so long as (f) ^ (j)Q at late cosmological times. The latter will 
introduce further restrictions on the allowed values of a, f3 and 5. With the exception of 
the k ^ —1 vacuum solutions, where there exists no FRW model within GR, we shall focus 
our attention on those particular parameter choices which reproduce Einstein's theory at 
late times and which we have just delineated in Eqs. ( |53[ ) - (|55|). When k = —1 the Milne 
model, a (X t, supplies the general relativistic solution at late times. 



V Theory 1: 2cj(0) + 3 



= 25i|l - 
constants. 



a > 1/2, Bi>0 



We study the evolution in the interval (p £ {0,(j)Q), allowing us to drop the modulus signs 
from 2uj + 3. Making the substitution u = (1 — (f)/(pQ), this bound becomes u £ (0, 1) and 
Eq. (|0D is 



^^^^d<p = Vm[ -V3lnK, a /I, 2, (56) 

(j) J M°/^(l — U) 

where K is an integration constant. In the models which asymptote to GR, u; ^ oo as 

A Vacuum solutions (A; = 0) 
1 Late-time behaviour 

For the theory to reduce to GR at late times we require u ^ as rj ^ oo. Choosing 
r/o = with A; = in vacuum we then have 

1 

-T, j. = ln{Krj), (57) 

A 1 — a/2 

as ti — > 0. A = \J?>l2B\ and its sign determines the sign of (2u;-|-3)^''^. The requirement that 
ti — > as r/ ^ oo demands we must have a > 2 and hence A > to ensure the right-hand 
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side of Eq. (^) is positive. Picking the appropriate right-hand side from Eq. (^) we thus 
obtain 



2 

2-a 







'r 









at large r/. This leads to the asymptotic form for t{rj), 

'a -2 



In2-Q {Kr])j 
In 2^ (Kf]) 



(58) 
(59) 



i(^)-^y^^'/'|l + ^ A 
The first-order inversion of this at large rj is 



2-a 2 



ln2-« [Ki]) 



(60) 



(61) 



This result is obtained by ignoring the factor in curly brackets on the right-hand side of 
Eq. (|60|). To obtain the next-order expression we substitute this first-order result into the 
weakly- varying ln(i^r/) term on the right-hand side of Eq. (|60|) , take the curly brackets onto 
the left-hand side and solve for ij. This yields 



vit) 



A 



a -2 



2-a 2 

ln2-« t 



(62) 



and 7/ ^ oo as t ^ oo. At late times we then obtain the solution 



a{t) 













CO 



2-a 2 

ln2-« t 



V20O/ \ 3 



a -2 



2-a 2 

ln2-« t 



(63) 
(64) 



2 Early-time behaviour 



At early times we see from Eqs. ( pOD and (56) that u ^ 1 as r/ — > is the only possibility 
consistent with u £ (0, 1) and A > 0. Treating n"/'^ ~ 1 in Eq. ( |5^ leads to the approximate 
relations 

ln(l-'u) = Aln(i<:7?), (65) 



and 



</)(r/) 



90 A 



K- 

— ■ 



(66) 
(67) 
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Taking the square root of the latter, integrating and inverting leads to 



\AK-^ 

When < A < 3 we see from the above that t — > as ry ^ and when A > 3 we have 
t — > — oo as ^ 0. In both cases we obtain the early-time behaviour 

3 2A 
A I 0r\ \ 3-A / 2 \ A-3 2A 



m ^ ) ,A-3 



'-A\ 3^ - ■ 



^ — ^ (70) 



When A = 3, Eq. (38) becomes 



A 

and t — > — oo as ?7 — > 0. The early-time evolution is 



Vit) = ew\\l^t\ , (71) 



0(t) ^ ^oK^exp|3W^t| , (72) 



a(.) ^ J^exp -^M!.,, (,3) 



as t — > — oo. Differentiating Eq. (|69D, we deduce that the early-time models will be expand- 
ing as long as A > 3 or A < 1. 

3 Minima 

We now probe the A; = vacuum models generated by this choice of coupling for the 
existence of expansion minima. Eq. ( p7| ) for these universes becomes 

a) • <^^> 

at a minimum. For u^: to be real for all a requires A > which is guaranteed when u G (0, 1). 
In general, however, we find 

< A < 1 ^ > 1 , (75) 
A = 1 ^ = 1 , (76) 
A > 1 ^ < 1 , (77) 

and when u G (0, 1), minima are only present when A > 1. We have the early-time solution 
for this model (Eqs. (^) and (^ )) in which a(0) = for A > 1 and the universe experi- 
ences a phase of contraction before bouncing and approaching late time general relativistic 
expansion. When A = 1, a(0) itself becomes a non-zero minimum (a^ — > A/K(j)o). 
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4 Exact solution 



The models with a = 1 are not described by the solutions presented above. These models 
tend (as (/> — > 0o) to the theory of gravity proposed by Barker [11|, for their evolution we 
find the exact results 



[K-^r]-^ + 1 



A 



(78) 
(79) 



When 7] ^ oo however, we observe that there is no combination of parameters allowing 
(j) ^ (pQ. Consequently, we shall not pursue this model any further. Theories with a = 2 



have been solved exactly and studied earlier by Barrow [12|, |52] and in ref. |45] 



When k 



B Vacuum solutions (k = —1) 
1 Late-time behaviour 

T the integral equation for the field evolution becomes (setting rjQ = 0) 

du 



at late times. As r/ — > co, tanhry 
may be approximated by 



W2(l -u) 

► 1 — 2e~'^^. Demanding u 



ln[Ktanhri] . (80) 
in this extreme, Eq. (^0|) 



1 ni-"/2 
Al - q/2 



IniT- 2exp(-2r/) 



(81) 



As ry ^ oo the right-hand side is finite and we require 1 — a/2 > (ie, a < 2) in order that 
the left-hand side does not diverge as n — > 0. We also require K = 1 to ensure that r] ^ oo 
and ti — > correspond to the same limit. The right-hand side approaches zero from below 
as 77 increases and u > 0; thus to keep the left-hand side negative we need A > 0. The 
form of y{r]) for the negatively curved vacuum models, selected from Eq. (p^), approaches 
Ae'^'^/A as ^ cxD and the late-time solutions are 



l-{X{2-a)r- 



exp 



a 



400 



1 + {A(2 - a)} 2-" exp 



4r/ 



a 



(82) 
(83) 



as 77 — > cx). Integrating and asymptotically inverting the second of these expressions we 
obtain the r]{t) relation 




{A(2-a)}^ 




(84) 
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Substituting this back into Eqs. (|2|) and ( |83[ ) yields the asymptotic forms 

l-{X{2-a)Y 



a{t) ^ 




(85) 
(86) 



as t — > oo. We note the asymptotic approach of this model to the general relativistic Milne 
universe at late times. 



2 Early-time behaviour 

At early times it — > 1 as ry — > and 

ln(l — li) ~ Aln [tanhr/] , 

which leads to 



(87) 



A 



-V 



l-A 



The functional form of this early-time behaviour is identical to that detailed in Eqs. ( |68| 
(|7^, after enforcing the choice K = 1. 

3 Minima 

Searching for stationary points in the scale-factor evolution for this theory, Eq. 
identifies 

cosh(2?7*) = An"/^ . (90) 

The cosh function is bounded below by unity, as is the value of when u G (0,1), 

which gives A > 1 as the condition for the existence of a stationary point. When < A < 1 
the universe expands monotonically away from an initial singularity at rj = . When 
1<A<3, a— >ooas7?— >0, the subsequent evolution will contain a minimum, allowing 
the initial contraction to 'bounce' back into late-time expansion. When A > 3 the evolution 
is expanding at all times. The special case a = 1 is treated by transforming r/ + % — > 
tanh(?7 -|- r/o) in Eqs. ( |7^) and (|79D. Solutions for the a = 2 case may be obtained from the 
solutions in ref.[B2| after applying the same transformation. 



C Vacuum solutions (k = +1) 
1 Late-time behaviour 

The positively curved {k = +1) vacuum solutions are governed by 

If du 

~XJ n°/2(l -n) 



In [K tan 7]] 



(91) 
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taking tjq = 0. In the GR hmit {u — > 0) this is approximated by 

1 ni-"/2 

In [K tan r/J . 



XI -a/2 

In conjunction with Eq. (|l^) we obtain 

'a -2 



(j){T]) 



1 



A 



2 

2-a 



ln2-« (Ktanr]) 



sin 2?] 



200 



1 



2 

{A(%^)}~ln^ [KtanT]) 



(92) 

(93) 
(94) 



These expressions approach GR when a < 2 as r] ^ tan~^(i^^^) and when a > 2 as 
rj — > mr/2 where n is an integer. The behaviour when a = 2 can be extracted from the 
k = treatment in ref. |52] by substituting 2(ry + ryo) — > tanh(?7 + ?7o). 

D Radiation solutions {k = 0) 
1 Late-time behaviour 

When A; = Eq. (22) leads to the asymptotic behaviour 



1 ni-"/2 



Ini^ 



2r?o 



(95) 



AI-q/2 r/ 

where we have picked 2Tr]Q = A to fix the arbitrary integration rjQ. The requirement 
that u — > as — > oo then demands K = 1 and a < 2. This upper bound on a, in 
conjunction with the lower limit implied by the requirement that u;'u!~^ — > 0, yields the 
powerful constraint 1/2 < a < 2. We examine models for which u £ (0,1), corresponding 
to i;^ G (Ojcjyo). In this case we have the constraint Ar/o > 0. Using y{rj) from Eq. (^), the 
late-time evolution of (p and a as functions of r/ is 



a)Ar/o 



V 



2 

2-a 



a [r]) — 
The latter of these allows us to deduce 



1 + 



(2 - a)Xr]o 
V 



2 

2-a 



vit) 



1/4 



.1/2 



1 /a 



4 Va- 1 



{A(2-Q)ryo}2- 



t)n \ 2{q-2) ^ 



and hence, the evolution as a function of cosmic time: 



ait) 



4r 



A\2(«-2)' 



,1/2 



1 + 



4(a - 1) 



4(/.o 



1 

2-a 



1 



(96) 
(97) 

(98) 

(99) 
(100) 
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as i — > oo, so there is a power-law approach to the GR solution at late times. 



2 Early-time behaviour 

At early times, we analyse the behaviour in the neighbourhood of u ^ 1, ie 



ln(l 



A In 



77 



V + 2% 



(101) 



The left-hand side of Eq. (101) tends to —00 as GR is approached. In order that the right- 
hand side approach the same limit we see that ry — > when A > and 77 — 2?7o when 
A < 0. Analysing the former case leads to 



(t>o 



-7] 



(102) 
(103) 



after setting 2Trio = A. The behaviour of this system as a function of cosmic time t mirrors 
that described in Eqs. ( |68| ) - (^) with replaced by (2ryo)^'^ and A replaced by 2rr7o. 
The bound Ar/o > implies r/o > when A > 0. When A < 0, — > — 2r/o at early times. 
The inequality Ar/o > implies % < and hence —2770 > 0. the evolution is approximately 



(104) 
(105) 



The t-parametrised behaviour of these equations is given by Eqs. ( pq ) - ([7^), after applying 
the transformations 



A 
A 



-2r/o)^ 
-A, 

-2rr?o . 



(106) 
(107) 
(108) 



3 Minima 

As with the vacuum solutions, we can search for turning points of when A > 0. 
Eq. (13) is 




(109) 



where subscript * denotes the value of a quantity at the stationary point. When < A < 1, 
the scale-factor expands away from a = at 77 = by Eq. (^). Monotonicity of r] implies 
r/* > and we know r/o > from the sign of A, which together give rise to the inequality 

u*/' > ^ . (110) 
A 

The range of < u < 1 together with the condition a > confirm A > 1 as a necessary 
condition for the existence of minima. When A = 1, r/=K = 0, and the universe expands 
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from a non-singular state of size a ~ {2rio)^^^T/(f)Q. When A > 1 the solution is initially 



contracting, bounces at ry* and tends to general relativistic behaviour at late times. Eq. ( 110 



) still applies when A < 0, and asserts that none of them can contain stationary points in the 
evolution of the scale- factor. This does not effect the models in which — 1 < A < 0, which 
begin from a = at r/ = — 2r/o and monotonically expand. When A < — 1 they are initially 
contracting and, due to the absence of minima, will always contract and never approach 
late-time general relativistic expansion. 



4 Exact solution 

The case a = 1 possesses a simple exact form, which is instructive. Solving Eqs. 
and (|2l|) leads to 



{r] + 2r]of + K^7]> 



2 ' 



Tr]{r] + r]o) (r] + 2r]of + K^rj^ 



(111) 



(112) 



where in these expressions we have again made the special choice 2rr7o = A. Examining 
the large rj limit reveals that (/) ^ (f'o only K = 1. In this case we integrate the 

asymptotic a(rj) to obtain 



and hence 



a{t) 



2 \A<poJ 



y^oJ t 



r 



valid as t ^ oo. At early times the solution becomes 

<^Wi - 

a\ri) - 



(113) 

(114) 
(115) 

(116) 
(117) 



as ?7 ^ 0. The behaviour of the i-parametrised version of this model is given by Eqs. ( |68| 



(73) after applying the transformations 



A 



(118) 
(119) 



and remembering K =\. Minima are a feature of this model when |A| > 1. 
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When k 



-1, Eq. 
1 



E Radiation solution (k = —1) 

61) becomes 



du 



InK = In 



(^2 _ p2^1/2g2(r,+r,o) _ p - ^ 



(^2 _ r2)l/2g2(»,+r,o) _ r + yl 



(120) 



1 Late-time behaviour 



Expanding the right-hand side at large i] and integrating on the left as u ^ leads to 
the approximate formula 

1 ui-°/2 2A 



-2rj 



121 



Al-a/2 (A2-r2)i/2 

where we have chosen A' = 1 to ensure ^ 0o at late times. The right-hand side of 
this expression tends to zero at late times and consistency on the left as u — > requires 
1 — a/2 > 0, or a < 2. If (pn > 0, as it must be when ^ c^o and < (p < (j)o, we recover 
the inequality XA > 0. Selecting the k = —1 form for y{ri) from Eq. ( [T^ ) we have the 
approximate solutions 



0(7/) 



AX{2 - a) 



(A2 - r2)i/2 



2 

2 — a _ Ar) 

e~ 2^ 



(A2_r2)l/2 
400 



,2r? 



. , r AX{2-a) \ 

l(A2 -r2)i/2 / 



2 

2 — a _ 4t? 



(122) 
(123) 



as 7/ ^ oo. Integrating and inverting Eq. (122) to next-to-lowest-order yields 
r?(t) 

y/±- — i - - 

AX{2 - a) 



(A2 -r2)V4 
1 f2 + a 



In 

"^2 \2-a 
t — oo as 77 — > oo and hence 

0(t) ^ 00 



2 

2-a 



(A2 -r2)V2 



rylA(2 - a) 



(^2 _ r2)l/4 



4 

" 2-a 



t ~~ 



(124) 



1 



400 



2 

2 — a 4 



a{t) 



t 



1 + 



2 f ^A(2-a) | — 



(2-a 



400 



(125) 
(126) 



as t — > cxD. Again, we observe power-law approach to the GR Milne solution at large t since 
a < 2. 
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2 Early-time behaviour 

In the early-time limit we find u ^ 1 and Eq. ( |120| ) becomes 

(^2 _ p2^)l/2g2(r,+r,o) _Y- A 



ln(l — li) — > Aln 



(^2 _ p2)l/2g2(»,+r,o) _ r + A 



(127) 



and as n ^ 1, ln(l — u) ^ —oo. To simplify the analysis we pick 770 such that n — > 1 as 
r/ — > , which requires 



=2% 



A + r\V2 



A-T 



(128) 



Combined with the k = —1 version of y(rj) from Eq. (^), this leads to the limiting forms 

'A + r 



(t>{v) 



A 



A /A + T 



A 



(129) 
(130) 



and the t-parametrised evolution will be that of the flat vacuum model described by Eqs. (^ 
- (73) after we apply the relabelling 



K 



A + T 



A 

3 Minima 



(131) 



To complete the study of this class of solutions, we search for points where the gradient 
of vanishes, Eq. (27) gives the condition 



^2 \ 1/2 



^2 



cosh [2(rj + 7]o)] ' 



(132) 



for a stationary point to exist. Since a > 0, u < 1 and cosh.[2(rj + 7]q)] > 1 this expression 
is equivalent to the inequalities 



.2 \ 1/2 



1--^ < A, A>0. 



I: 



1/2 



> A, A < 0. 



(133) 
(134) 



Since A"^ > E^, we know that < (1 — r2/yl2)i/2 ^ Hence, models with A > 1, in which 
a begins collapsing from infinity, will bounce to mimic GR expansion at late times. 
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F Radiation solutions {k = +1) 



If we selecting the necessary right-hand side from Eq. ( p^ ) for the finite, closed-universe 
models with k = +1, then Eq. (^) is 



du 



Xj W2(l-ti) 



In 



K 



'rtan(r? + r?o) + {T'^ + A^^/"^ - 
r tan(r? + r/o) + (r2 + ^2)1/2 + ^ 



(135) 



As u ^ 0, we obtain the following approximate expressions for the behaviour of 4> and a: 



(pii]) 



which are valid as 



2 

2-Q 



X ln2-a 



1 

200 



K 



X ln2 



'rtan(r? + ??o) + (r2 + ^2)1/2 
rtan(7? + 7?o) + (r2 + ^2)1/2 + 

|r + (r2 + ^2)i/2g^2(r/ + 7?o)]}|i 

rtan(r? + 7?o) + (r2+yL2)i/2 
rtan(?7 + ?7o) + (r2 + ^2)1/2 + ^ 







(- 


H¥)l 







2 

2-Q 



(136) 



(137) 



G Perfect-fluid solutions {k = 0) 



We shall analyse the late-time behaviour of Theory 1, specified in Section IV, by the 
device of using the solution defined by the choice of field evolution 



(/>(0 = </'oexp(^ln^^) , 
where A and B are constants. This gives rise to the generating function 



AB 



(138) 
(139) 

From Eq. ( |14C| ) it is clear that to keep the left-hand side positive we require the combination 
(C - 1)^-'^AB > 0. This leads to 

e 



and the scale-factor 



,3(2-7) 



3(2-7) 
^0 



e2ln^-^eexp(^ln^e) . 



and hence 



2lj{C) + 3 



4-37 



4-37 e 
2 



3(2-7)^ 



^ In^-^ £ + ^£2 



(141) 



(142) 



We require </) — > 00 at late times and this occurs as £ ^ 1 when B > and as ^ cxd when 
B < 0. When B = 0, (p = (jioe^ and when A = 0, <p = <po; in both cases the theory is GR 
at all times. 
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1 Late-time behaviour 

When < B < 1, we find from Eq. (|2|) 

l-B 

t ocln ^ , 



(143) 



as ^ — > 1. The power of the log is always positive for < 7 < 2 and t — > in this limit. Since 
we are only concerned with limits corresponding to large t we exclude the range < -B < 1. 
When B > 1 the coupling function tends to 



2co(e) + 3 - ^-^^ ln-(^+i) £ 

^^^14) +3^ 3^2-^)2 AB ^' 



(144) 



and Eq. 



becomes 



dt 



l-B 



3(7-1) 
^0 



4-37 

3(2-7) iAB<p,f^ V^45" 
As £ — > 1 this tends to 



1/2 



£^U-7i In 2(2-7) £exp 



7-1 
2-7 



In^C 



1 - s a; 



3(7-1) 




3(2-7)(AS0o)^ ^ 



4-37 



1/2 37-4-B7 

(£ - 1) 2(2-7) , 



and 



where 



7(1-B) 



37 (A5(/)| 



7-1 
I 2-7 



a. 



3(7-1) 




4-37/ 



Eventually, this leads to 

m 
m 

ait) 



2(2-7) 2(2-7) 

1 _)_ jp-,(i-B) tia-B) ^ 

2B(2 — t) 2B(2-7) 
1 + AF 7(1-S) t 7(1~S) 



i?2/37 



.^2/37 



iAB<Po 



3(2-7) 



2(2-7) 
2F'y(l-B) 2(2-7) 

1 H ; ^t^(T^ 



3(2-7) 



(145) 
(146) 

(147) 
(148) 



(149) 
(150) 

(151) 



as t — > 00, and there is power-law approach to the GR solutions in this limit. The coupling, 
in terms of the field, approaches 



2u;{(j)) + 3 



3(2 - 7)2 B 



00 



(152) 



The negativity of d.^ / dt implied by Eq. ( |145| ) for B > 1 implies that £ approaches unity from 
above. As we noted earlier (after Eq. (|140|)), when £ > 1 positivity of B implies positivity 
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of A and thus Eq. (138) confirms that (/> — > 0o from above in these theories, ie i;^ G {(po, oo). 
When B < the couphng function, as ^ — > cxd, is 



2u;(0 + 3 



and the temporal line element becomes 



dt 



3(7-1) 
^0 



3(2 - 7)2 AB 



4_3^\l/2 2(7-1) j(l-B) 

* ^2-7 In 2(2-7) ^d^. 



(153) 



3(2 - 7) ^ABc^o 



2-7 ^ 



Integrating, this gives 



where 



7 tjl-B) 

i(0 - C~^C~ln^(^ ^, 



C 



37 {ABcp, 



7-1 
I 2-7 



a, 



3(7-1) 




4-37 



Inverting asymptotically in t, this gives 

2-7 2-7 B-l 2-7 

£^{t) t '< In 2 t 7 J 

and hence there is logarithmic approach to the GR perfect-fluid solutions 



m 

ait) 



1 + ^ln^t 7 



2-7 



ao 



C2/37 



_^2/37 



1 + — n^ t 7 

3(2-7) 



as t — > CXD. The late-time behaviour of the coupling as a function of (j) is given by 

4-37 4A^i/B i_s 

2uj(4>) +3^ — '-^ In— 

3 2 - 7)2 B 



(154) 

(155) 
(156) 

(157) 

(158) 
(159) 



(160) 



At large t, ^ ^ oo and must do so from below to maintain positivity. Consequentially, 
B < implies A < and from Eq. (|138|) cp ^ (po from below, ie (p G (0, (po) in these models. 



2 Early-time behaviour 

When B > 1, a(^) has a zero as ^ ^ iff {—1)^A < 0. At late times ^ — > 1 from 
above; however, this is a manifestation of the fact that (,{t), for the exact theory defined 
by Eq. ( |138| ), is not monotonic. That this is correct may be demonstrated using Eq. ( ^2|) 
and solving \/2cj -|- 3 = for ^ when B = 2 then further showing that \/2ci7+~3^ 7^ at 
that point. As ^ ^ 0, In^ < and {-l)^-^AB > so that a{C) > 0. When B > 0, 
{-1)^-^A > and hence (-1)^^ < 0, ie the universe is always singular at ^ = in i3 > 1 
theories. In this limit 

(4-37)2 



2uj{0 + 3 



3(2-7)2 



(161) 
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ie a BD theory. 

When < 0, a(^) has a zero as ^ — > 1 iff ^ < 0. This arises since \/2uj + 3 > as ^ 1 
and ^ approaches unity from above (In^ > 0). As in the preceding paragraph, a(^) > 
imphes AB > 0, B < enforces ^ < and the evolution begins from a singularity at ^ = 0. 



Again, the form of the coupling function at early times is given by Eq. (161), BD theory. 
We will not present explicit solutions showing the approach to BD theory since these 



models are not early-time limits of the theories, defined in Section |V|, that we are interested 
in; they are merely early-time limits of other theories which happen to asymptote (at late 
times) to the theories with which we are concerned. 

We now highlight some special cases of these models for particular values of 7. 



3 Dust models 

These arise by substituting the choice 7 = 1 into the asymptotic relations already 
derived. When i? > 1 we find 



ait) 



3\ l-S 1 B 2B 

A^-B B^-BU-B 



^\^^^ "0 f2/3 



l+B 

l+(-j {AB)—t— 



When i? < the solutions tend to the GR solution only logarithmically 

a{t) ^ 

All of these expressions are valid as t ^ 00. 



l + A\n^t 
Qq ,2/3 

2 J ^;/'' 



(162) 
(163) 

(164) 
(165) 



4 Inflationary models 

Inflationary models driven by a false vacuum equation of state may be derived from the 
choice 7 = 0. Although there are varieties of inflationary universe with — 1/3>7>0, and 
these can easily be found from the formula for the general 7 solutions given above, we shall 
confine our attention to the 7 = case which is not described by the previous formulae. 
It offers an excellent approximation to many slowly changing scalar-field potentials. In 
this case we can view the scalar-tensor coupling as providing a second scalar field, thereby 
offering the chance for double infiation to occur. However, it is not sufficient simply to 
substitute 7 = into the existing expressions since the qualitative structure of the solutions 



is different. When B > \, the form of 2a; -|- 3 is as presented in Eq. ( 144 ) and the temporal 
line element of Eq. ( |145|) can be approximated by 



dt ~ 



{l-B) 



3a3 



.1/2 

-^In-i^dlne, 



(166) 
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as ^ — > 1. Integrating this expression we find 

which — > oo as ^ ^ 1. This can be inverted, yielding 

^{t) ~ exp < exp 



(167) 



(168) 



and hence, 



a{t) 



1 + A exp 



3alBt 



(1 - B)<P'J' 



(^i?0o)'/''''^V2('Ao)i/2 



1 f 3ah 
1 + :j exp <^ -j^ 



(169) 
(170) 



as f ^ oo. Here we see exphcitly the possibiHty of double inflation arising from the sequen- 
tial effects of the (p field and the p = —p stress. If -B < 0, t ^ oo as ^ ^ oo and 2(a; + 3 is 
given by Eq. ([1531). The differential S, — t relation, Eq. (|154|), is then well approximated by 



dt 



200' 



1/2 



3a^ V 2 



A 



exp I -— In^n dln^ 



(171) 



Making the substitution ( = In^ ^ we can integrate the above equation to obtain 



1/2 



200 
3a^ 



■In^ 



1 



A 



2(B + r 



■ln«^ 



(172) 



as ^, and hence t, tend to infinity. Asymptotically, we obtain 
^(t) exp 



3a|]t ^ 

—[/2 



m 

ait) 



Mo 

l + ^^^"o 



200 ' 



ao 



3a; 



l-B 
3 \ ~6~ 

\ 



[ABc^of' \2ctl'\ 



t 6 exp 



' apt ^ 



(173) 
(174) 

(175) 



as t ^ oo. 



5 The Connection to the parameters of Theory 1 

We now use the results derived earlier in this section to model the late-time behaviour 
of Theory 1. Consider first universes in which (p ^ (pQ from above, ie G (0o, : oo). The 
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approach of Theory 1 to the relativistic hmit in this direction can be accurately modehed 



using the theory defined by Eq. ( 138 ) with A > 0, B > 1. From Eq. ( [15 



2a;((/>) +3 



4-37 (l-i3)^ i/B 
3(2-7)2 B 



— -1 



B J 



(176) 
(177) 



as (/> — > (/>o. This expression is essentially the same as definition of the Theory 1 coupling for 
(j) £ (00) oo) introduced in section 4. Explicitly, we may obtain the asymptotic behaviour of 
Theory 1 by making the following identifications between its parameters and the parameters 



of Eq. ( pSl) : 



A 
B 



6^1 (2-7)^(a-l) 

(4-37) («-2)' 
1 



1 

a-l 



a 



1 



(178) 
(179) 



The constraint on the S > models such that they approach GR at late times, namely 
> 1, is equivalent to —{B + 1)/-B > —2. This is a very restrictive condition because the 
function —{B-{- 1)/B is naturally bounded above by the value —1. Thus, for perfect-fluid 
universes with < 7 < 4/3 Theory 1 can only be expected to converge to the general 
relativistic value of (f) from above if 1 < a < 2. When (j) converges to (pQ from below, ie 
<j) G (0, 0o) we can approximate the late-time behaviour of Theory 1 using the solutions for 
A < 0, B < 0. The asymptotic form of the coupling in this case is given by Eq. ( |160D 



2a; ( 



+ 3 



4-37 4(-A)-i/^ 
3(2-7)2 (-B) 



1-^ 



l-B 
B 



(180) 



as (f) ^ (1)0. The behaviour of Theory 1 at late times may be found by substituting the 
expressions 



A 
B 



SB, (2-7)^(-l)' 
2 (4 - 37) (1 - a) 



1 — a ' 



(181) 

(182) 



into the formulae describing the asymptotic evolution of the theory deflned by Eq. (138). 
The choice B < leads to the constraint < — 1, which becomes a > 1 for Theory 1 as 
d) ^ <j)n from below. 



VI Theory 2: 2uj{(J)) 



3 = B2|ln(0/0o)| 
constants. 



-25 



; S > 1/4, ^2 > 



The left-hand side of Eq. (20) for this choice of the coupling function, when cp £ (0, (po), 
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IS 



/ 



(2u;(0) + 3) 



1/2 



-dcp 



1-5 



/Kin 



(183) 



where K is an integration constant. We now investigate vacuum and radiation solutions for 
this case, under the assumption that (/> G (0, 0o)- 

A Vacuum solutions (A; = 0) 



Selecting the zero-curvature right-hand side from Eq. (|2[)D and the necessary form for 
y{r]) from Eq. (0), Eq.( |183| ) for 5^1 reads 



ln( ^ 



l-<5 



-A2(l-<5)ln[i^(r? + ?7o)] , 



(184) 



where A2 = \/?>/B2. The left-hand side is positive, for the right-hand side to follow suit as 
r/ — > 00 requires A2(l — 5) < 0. As r/ — > 00, the right-hand side^ 00 and we must have 5 > 1 
for this to occur on the left as (j) ^ This implies A2 > and we obtain for the field and 
the scale-factor, when 5 7^ 1, the exact expressions 



<^(r7) = (/.oexp[-{A2(<5-l)ln[if(r/ + r/o)]}~J , (185) 
a'(^) = ^(r/ + r?o)exp{[A2(<5-l)ln[K(77 + 7?o)]]T^} . (186) 

When 5 = 1, the conformal time-evolution of the field and the scale- factor are given by 



0(r?) 



lexp [-K{r] + iqQ)] 



-A2 



— (?? + m) exp ] [K{ri + riQ)Y ^ \ 



1 Late-time behaviour 



(187) 
(188) 



As ?? ^ 00 the late-time behaviour can be modelled by Eqs. ( ]60| ) - ( p^ ) under the 
simultaneous transformations: A ^ A2, a — > 25. Setting r/o = for simplicity, we see that 
the scale- factor tends to zero as ry — > when 5 > 1. When 5 ^ 1 we require A2 > if (/> is 
to tend to (/>o as — > 00. We find 



(Pit) ^ (Po 



l-K 



-\2 



AA 



-A2/3 



-2A2/3 



(189) 



which approaches (Jjq as t — > 00. The asymptotic form of the scale-factor is 



a{t) 



/3£y/3 



/3 



-A2 



V20O/ 

with limiting behaviour a oc 

^1/3 

as t ^ 00. 



1- A2 
3-2A5 



AA 



-A2/3 



-2A2/3 



(190) 
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2 Early-time behaviour 



Examining the form of a(ry) when 5 ^ 1 in Eq. ( 185^ , we see that the choice % = in 
ensures a(0) = 0. As rj ^ 0, \nr] ^ —oo and the exponential factor in a either tends to zero 
or a constant, depending upon the sign of 1 — 6 . However, the bound 1 — 6 < on the power 
of the logarithm in the exponential guarantees that we always obtain a(ry) oc \fAj^r^l'^ 
as r] ^ and hence r? (9(/>o/4^)^/^t^/^. This final expression implies t ^ as 77 ^ 
and the early-time formulae for (j){t) and a{t) are identical to Eqs. (|60|) and (|6^) , under the 
transformations A ^ A2, a ^ 26. 

When 6 = 1, the dominant behaviour in a{r]) as ?? — > is contained in the exponential, 
from which we may conclude (after setting r/o = 0) that 

1 



{KrjY 



(191) 



This may be inverted asymptotically in ry^^. To first-order we obtain rj ~ K^^(2\iit)^^^'^^ 
as r] and t tend to zero (noting the sign of A2 and the necessary monotonicity of rj{t)). The 
next-order corrections at early-times follow by substituting this lowest-order result into the 



weakest dependence in Eq. (187 ), ie the power-law factor, and solving for rj. This yields 

-1/A2 



T] 



K 



-1 



lnt+ 1 + 



3 
2X2 



In(lnt) 



and hence 



m 

ait) 



"2 In ^2 ^t. 



oc 



t (In t) 



1+1/A2 



(192) 

(193) 
(194) 



as t — > 0. 



3 Minima 

Examining the 6^1 vacuum models for the presence of minima we find 



■exp{-[A2((5-l)ln(i^7?)]T^} 



(195) 



For the exponential to tend to zero we require its argument to tend to —00. the bound on 
6, namely 6 > 1, implies that rj ^ K for this to happen. At this point {Krf) 
00; nevertheless, this logarithmic divergence will be insufficient to counter the exponential 
convergence of the pre-factor. We can also expect to see turning points in the evolution of 
a when the factor in square brackets vanishes. This happens when r? = jy* , where 



K- 



exp 



(-A2)-^/^ 



1 



(196) 



One can show that the conditions for the factor in square brackets in Eq. ( |195| ) to vanish 
also guarantee that the exponential will be well behaved. When 5 = 1 we have 

A 



(a )^ 



1 + A2-fC~^'ry-^2 



exp 



(197) 
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Krj > and so the exponential will vanish when rj ^ 0. When this happens the factor in 
square brackets diverges, however its divergence is quashed by the rapid convergence of the 
exponential, confirming the presence of a stationary point at r/ = 0. There will also exist 
a stationary point in the evolution of a when the square-bracketed factor itself vanishes in 



Eq. (197). This occurs at r/*, where 

V* = i-X2y^^'K-\ (198) 
At this point —{Kvi)~^'^ = A2 ^ > and the exponential is well behaved. 

B Vacuum solutions {k = —1) 

Similarly, we analyse the behaviour of the negatively curved models. Selecting the 
k = —1 versions of Eqs. ( |20| ) and (0), we obtain for the field and the metric, when 6^1, 
the expressions 

(l){r]) = (/.oexp{-[A2(5- l)]^lnT^ [Ktanh(r? + ?7o)]} , (199) 

"^(^) = TT- sinh [2 (r/ + r/o)] 
200 

xexp| [A2((5- 1)]^^* InT^s [Ktanh(r? + r/o)] }. (200) 

This solution will approach the Milne model at late times under the conditions < S < 1, 
K = 1. The approach of (p to (po from below requires A2 < 0. When 6 = 1, setting 770 = 0, 
we have 

(f){rj) = (/)oexp{-K-^nanh-^2^| , (201) 

c? (r/) = — sinh(27?) exp <^ K"^^ tanh'^^ ^ I _ (202) 
200 1^ > 

Because there is no choice of K which allows to tend to 0o at late times we do not pursue 
this model any further. 

1 Late-time behaviour 



The late-time behaviour is given by Eqs. (|8J) - ( pSD under the substitutions A = A2, 
a = 25. 

2 Early-time behaviour 

Picking r/o = 0, = 1 in Eqs. (199) and ( |200| ), we have as r/ ^ 0, 



0(r?) ^ 0oexp{-[A2((^- l)]^lnT^??} , (203) 
«^(^) 4-^exp|[A2((5- l)]T^ln^r/| . (204) 
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We know < 6 < 1 and hence 1/(1 — 6)>1 and a^irj) will be dominated by the exponential 
factor as r/ ^ 0, and hence a ^ 0. Integrating and inverting a(r/), we obtain 



A,(^-l) 



l+3(^)ln-H) 



(205) 



where t — oo as rj ^ 0. Using this relation the early-time behaviour is 



<P{t) ^ </,of-2exp|-3(^^-^^^lni-^(-t)| , (206) 
2.x f 2i-^lni-^(-t)l , , 

" ^ fc^-H A.(i-1) ' (^0^) 



as t — > — oo. 
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Differentiating Eq. ( pOO] ) yields 

(a^) = 4- (cosh(27/) + (-A2)i^ (1 - 5)1^ lnT^(tanhr?)| 

xexp|[A2(l-(5)]^lnTi5 [A"tanh(?? + ??o)]} • (208) 

The exponential cannot tend to zero, since its argument is always positive. If we examine 
the pre- factor we find stationary points exist at 7/ = r/*, where 

cosh(2r/*) = -{-X2)^{l-S)^ [In (tanhr?^)]^^ . (209) 

This is not soluble analytically, although we may gain a bound on its value by demanding 
that the cosh function always be greater than unity. We obtain 

r?* < arctanh jexp ^- ^^ | . (210) 

C Vacuum solutions (k = +1) 

When (p lies in the range < (p < (f)Q, and 6 ^ 1 we have the exact solution 

cl){r]) = 0oexp{-[A2(5-l)]^lnT^ [Ktan(?? + ??o)]} , (211) 
a\r]) = ^sin[2(r/ + r/o)]exp{[A2(5-l)]^lnT^ [Ktan(77 + 77o)]} . (212) 

In the particular case 6 = 1, we have instead 

(Pirj) = (t>oexp[-K-^Han-^^{i] + i]o)} , (213) 

a^(^) = 77— sin[2(r? + ryo)]expi-i^-^nan-^2(^ + j^Q)l ^2U) 
200 
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D Radiation solutions {k = 0) 



Fixing the origin of conformal time in Eqs. (22) and ( |2lD such that 2Tr]Q = A, we obtain 
the exact results, for 5 / 1 



(pi-q) = (/)oexp <^ - [A2(5 - Im-^ 



K 



7] 



a\r]) = — 7?(r/ + 2r/o)exp|[A2((5- hii 
In the special case of (5 = 1, we have the exact relations 



T] + 2?7o 

1 1 

K 



7] 



7] + 27]o 



<p{7]) 



90 exp 

r 





Kri 






7] + 27]o 





»0 



r/(ry + 2r]o) exp 



Kr] 



7] + 2?7o 



(215) 
(216) 

(217) 
(218) 



To ensure > 0, we require r]{7] + 2r/o) > and the modulus signs in the above expressions 
can be dropped. 



1 Late-time behaviour 

The asymptotic behaviour may be obtained from Eqs. ( [9^ ) - ([99|), by applying the 
transformations a — > 25, A ^ A2. When S 1 Eq. (|215| ) serves to bound the allowed 
parameter values in order that (j) ^ (po ss 7] ^ 00. We can see that we require 6 < 1 and 
K = 1, delimiting the allowed range of(5tol/4<(5<l. We take A2% > to ensure both 



sides of Eq. (22 ) exist on the correct domains. Examining Eq. (217), we see that when 
6 = 1 there is no non-trivial choice of K which permits — > 0o at late times. Therefore we 
exclude it from further analysis. 



2 Early-time behaviour 

At early times when 5^1 the behaviour is harder to ascertain and we need to make 
use of logarithmic approximations. From Eq. (l216| ), by demanding that (/> G (0, i;^o)> we see 
that the lower limit of r] occurs as 



(i) 7/ ^ , A2 > 0, ?7o > 0, or , 

(ii) 7] -27]Q , A2 < 0, r/o < . 



Case (i) gives 



a(r/) ~ 7] 



i/2exp|i[A2(5-l)]^lnT^ 



(219) 
(220) 
(221) 

(222) 



which, for the range of values of 5 to which we are confined, is dominated by the exponential 
at small 7]. Integrating this expression and inverting approximately then leads to 



7]{t) = 2r/o exp 



A2(l - 6) 



In(-t) + 3 



A2(l-<5) 



(223) 
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and so i — > — oo as — > 0. The field and the metric are then given by 

(j){t) ~ t-^expjs ^ '^l \. ln^-^(-t)| , 



a{t) 



t exp 



A2(<5-1) 

2i-5 



as t — 
r] ^r] + 2r]o, A2 



A2(l-5) 

00. Case (ii) can be modelled by applying the transformations: 770 



(224) 
(225) 



-A2, in this order, to Eqs. ( p23|) - (|224| )- 



3 Minima 



Since the scale-factor is infinite at early-times, we know there must exist at least one 
minimum in its evolution, in order that we obtain general relativistic expansion at late 
times. Differentiating Eq. ( |216| ), we find 



2r 



V + Vo + Vo (-A2) i--* (1 - 5) Ini- 



7] + 2r]o 



X exp <^ [X2{S - l)]!-* Im- 



, , - (226) 

For the parameter choices we are confined to, the exponential factor in the above expression 
is a monotonic function, existing in the range (1, cxd). Thus we search for zeros of the pre- 
f actor, finding them to exist at rj^:, where 



??*+??o + ??o(-A2) 1-^(1 - 5)^-s Ini- 



V* 



0. 



which is non-analytic and must be solved numerically for particular 6, A2 and r]Q . 

E Radiation solutions (k = —1) 

For the negatively-curved models with 5 7^ 1 we have exactly 



(227) 



(j){r]) 



200 



>o exp 
1 

X Ini-* 
1 



1 

1-5 



A2(<5-1 

r tanh(?? + 7?o) + (r^ + AY^^ - A 



K 



K 



xlni- 

When 5 = 1 these become 

(t){ri) = (/>oexp<! 



r tanh(7? + T/o) + (r2 + ^2)1/2 + A 
T + {T"^ + A^fl^ sinh [2(7? + 7?o)]] expj- [A2(5-1)]t' 
rtanh(?? + ??o) + (r^ + A^y/"^ - A 



(228) 



rtanh(77 + 770) + (r2 + ^2)1/2 + A 



T tanh(r/ + ryo) + (r^ + ^2)1/2 _ ^ 



rtanh(r? + r?o) + (r^ + ^2) 1/2 + ^ 



(229) 



(230) 
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J_|_r+ M2 _ p2Nl/2g2(r,+r,o)| 



exp J K~^^ 



rtanh(r/ + ryo) + (F^ + ^2)1/2 _ ^ 



rtanh(r/ + r/o) + (r^ + ^2)1/2 _^ ^ 



(231) 



1 Late-time behaviour 



The late-time behaviour may be derived from the Theory 1 solutions, by the simulta- 
neous transformations A — > A2, a — > 25. When S ^ 1, Eqs. (124) - ( |125| ) require AX2 > 
and 5 < 1 to ensure late-time approach to GR. When 6 = 1, there exists no value of K that 
will admit a tendency of (p to (po at large r], and for this reason we pursue these models no 
further. 

2 Early-time behaviour 

We make the simplifying choice for the origin of r/-time, 



=2»70 _ 

\A-T 

Eq. ( |228| ) for the 6^1 solution then becomes 



(232) 



Hv) - <l>o exp 



A2(5-l)ln 



r + ^ 

A 



1 



(233) 



as T/ ^ 0. When A2 > 0, — > as ry — > 0. When A2 < however, there only exist solutions 

1 

when (—1) is real. For the limiting form of the solution, we find 



r,{t) oc t^/\ 

(p{t) oc exp I 

a{t) oc t^/^ , 
as r], and hence t, tend to zero. 



^(5-l)lnt 



(234) 
(235) 
(236) 



3 Minima 



We define 



such that 



e^^-1 

=2»? _ p-4?7o ' 



(237) 



^|(^2_r2)i/2^^g^P(r? + r?o)] 
+ [-r + {A^ - r2)V2 sinh [2(ry + %)] 
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(238) 



At a stationary point in the evolution of the scale-factor we require the expression on the 
right-hand side of Eq. ( |238D to vanish. Since a zero of the exponential would require (f) oo, 
which is outside our range of consideration, we search for zeros of the pre-factor. Minima 
exist at r]^:, given by the implicit formula 



{A^ - r2)i/2 cosh [2{r], + r/o)] + \-T + {A^ - T^)^/^ sinh [2(r/, + %)]' 



X (-A2)i-* (1 - (5)1-5 



22^* (1 



(e2r,. 



-4»7o^2 



i5 



=2r,. 



p2r?. 



-4r?o 



F Radiation models (k = +1) 

For the closed models with S ^ 1 we have 



lexp 



X lni-'5 



K 



[X2{S-1)] — 
r tan(r/ + r/o) + (T^ + ^2)1/2 _ ^ 



a'(r/) 



X In I-* 

and when 6 = 1 these become 



r tan(r/ + r/o) + (r2 + ^2)1/2 + ^ 
_L (r + + ^2)1/2 ,i^[2(r? + %)]) exp I [A2(5 - 1)]T^ 
r tan(r? + rjo) + (T^ + ^2)1/2 _ A 



K 



(j){r]) 



I exp 



X In 



K 



r tan(r7 + %) + (r^ + ^2)1/2 + ^ 
exp ^— A2 

r tan(r/ + r/o) + (r^ + A^/'^ - A 



r tan(r/ r/o) + {T^ + vl2)i/2 + A 
_L (r + + ^2)1/2 gij^[2(r/ + %)]) exp jexp 
r tan(r/ + %) + (r2 + ^2)1/2 _ ^ 



X In 



r tan(r/ + r/o) + (r^ + ^2)1/2 + a 



G Perfect-fluid solutions (k = 0) 



. (239) 



(240) 



(241) 



(242) 



(243) 



When 



bo from above, ie (j) G 



00) the approach of Theory 2 to the general 



relativistic limit can be accurately modelled using the theory defined by Eq. ( 138 ) with 
A>0, B>1. Recall Eq. (Q, 

\2 



2uj{(f>) + 3 



4-3-, 

3(2-7)2 B 



(244) 
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as (j) ~^ 4'o- Thus we can deduce the late-time evolution of Theory 2 by enforcing the 
relations on the perfect-fluid solutions presented in Section ^ 



A 
B 



3B2 (2 - {26 - 1) 



4(4-37) (s-i) 
1 



1 

2(5-1 



26 - 1 



(245) 
(246) 



when (j) ^ (pQ from above. Theory 2 can only tend to GR at late times with (f) going to (f)Q 
from above ifl/2<5<l, as follows from the constraint B > 0. When (f) ^ (po from below 
the behaviour of the coupling is given by Eq. (|16C1|), 



2uj((h) + 3 



4-37 4A"i/B i-a ^ 
In s — 



3(2-7)^ 



B 



The late-time behaviour of Theory 2 when (j) is in the range (j) £ (0, 
Eqs. (138) - ( |173| ) under the substitutions 



A 
B 



3(2-7)^^2 
4 (4 -37) (1-25) 

1 



1 

2S-1 



1-26 



(247) 
(248) 

is then given by 

(249) 
(250) 



The bound B < leads to the constraint 6 > 1/2, hence there exists a wide spectrum of 
models with < 7 < 4/3 in which (p ^ (pQ from below. 

VII Theory 3: 2u{(J)) + 3 = ^3 |l - ((/)/0o)^|"\ /3 > 0, ^3 > 

constants. 



For this choice of the coupling function, we find 



{2uj + 3) 



1/2 



In 



1 - y/U2 



where < 6 < (po and 



U2 = I 



A Vacuum solutions (k = 0) 

Starting with the flat k = models with G (0, 0o), we obtain 

1 - ^/u^ 



A2/3 



(251) 



(252) 



(253) 
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with A2 = y/S/B^ and fixing 770 = 0. Using Eq. ([l9|) we can then deduce the evolution of 
the field and the metric to be 



4^/1^ (bo {Krif^ 



X2I3 



1 + (Krj) 



2/(3 



1 + (Kr]) 



A2/3 



4l//3(/,o 

1 Late-time behaviour 



2/(3 



(254) 
(255) 



Examining the form of Eq. ( p54| ) we see that there is no choice of K for which 
(j) ^constant at large rj, precluding any possible approach to GR at late times. In spite of 
this, we find for the asymptotic behaviour 



(pirj) 



At late times we find 

vit) 

m 

a{t) 



1/2- 



\AK\ 



A2I 



3+1A2I 



3+1A2 



/3+ A2 



2AV2 



^3+|A2l 



2|^2l 
"3+|A| 2IA2I 



/AK|A2|\ 3+[A2i /3 + lAalxirlifT 1+1^21 



(256) 
(257) 

(258) 

(259) 
(260) 



Since there is no late-time approach to GR, we do not pursue the early-time behaviour or 
probe for the existence of minima in these models. 

B Vacuum solutions {k = —1) 

Similarly, for the k = —1 cases we have 



1 - Vu2 
1 + 



i^^2/3tanh^2/3 



V: 



(261) 



and convergence to GR as r/ ^ 00 demands K = 1. By direct comparison with the solution 
for fiat models we obtain the exact results for the evolution of 6 and a 



tanh^2 



V 



(1 + tanh^2/3j^j^/^ 

A (1 + tanh^2/3^^^/^ 

— smh{2r])- v 

2~6n tanh ^ r] 



(262) 



(263) 
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1 Late-time behaviour 

Asymptotically, these relationships tend to the pair 

"l - \lPe~'^ 
1 I A 



4>{v) 



1 + 



(264) 
(265) 



as r] ^ 00, and we find 



r]{t) ~ In 



leading to 



ait) 



+ 2402 



(266) 



(267) 
(268) 



as t — > 00, and so there is power-law approach to the Milne universe of GR at late times. 

2 Early-time behaviour 

Setting r/o = 0, we can expand Eqs. ( p62D and ( |263 ) about rj = 0. We find 

(j){r]) ~ 4^/^(/.or/l^2| ^ 



a2(ry) ~ 4-1/^-r/i-l^^l 



The latter of these leads to the early-time r](t) relation, 



2(1-/3) / (An \ 3-IA2I 2 2 

rj{t) = (r^) ' (3-|A2|)^^t^^ 



(269) 
(270) 



(271) 



Studying the form of the exponent in Eq. (271) reveals that when IA2I < 3, t ^ as 7/ — > 
and when IA2I >3,t— >— cxDas?7^0. The early-time evolution of the universe as a 
function of t is then given by 



2(3-|A2|/3) 



2IX2I 2IA2I 



0(i) ^ 2 '3(3-1^2!) ,/,q"-i^2I ^lAiF3 (3 _ |A2|)^^^ t^=^ , (272) 

/3(|A2|-l)-2 / 4 \ TgrjATT 1-I^2l 1-I^2l 

^ 2 /3(3-|A2i) _ (3- |A2|)^^t^^ , (273) 

\(Po/ 

as t ^ 0, —00 accordingly. When IA2I =3, the t -dependent evolution is given by 



r?(t) 



exp 



(274) 
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exp 



A 



ait) 



exp 



(275) 
(276) 



as t ^ — oo. 



3 Minima 

Differentiating Eq. ( p63D , we find 

(277) 



y times, ie 

implies 



a2 



2A 

cosh(2ry) H ^-g A2 



sinh(2r7) 

At early times, ie as t/ — > 0, ^ if IA2I < 1- Since sinh(2?7) ~ 2r/ as ^ 0, Eq. ( 269) 



oc ri-\^^\ [. . .] , (278) 



as r/ — > 0, and thus the power-law pre-factor in (a^)^ can never vanish at early times. 



The field (f> evolves monotonically, by Eq. ([T7|), and tends to a constant at late times, which 
guarantees that (prj cannot diverge during the subsequent evolution. Eq. ([T^), in conjunction 
with the bound < (p < (f)o, then ensures that there can be no further zeros of a as the 
universe evolves. Eq. ( |264| ) ensures (a^)^ is non-zero as ^ 00. In general, stationary 
points in the evolution of a arise at r/,,, given by 

2A 



obtained from the square-bracketed factor in Eq. ( ^771 ). Remembering the positivity of both 
cosh(2ry^) — 1 and 77*, we can derive the interesting result that when IA2I > 1 the location 
of the minimum is constrained by 

T]^ < arctanh , (280) 



.A2 + I 

and when IA2I < 1 stationary points do not exist. 

C Vacuum solutions (k = +1) 

Finally, we note the existence of k = +1 closed-universe solutions, governed by 

^-^ =K^2f^t^n^2Prj. (281) 

1 + VU2 

This leads to the exact solution, in conformal time. 



4^/^(t>oK^^ tan^2 



(1 + /s:A2/3tan^2/3^) 



2/13 



js^ fl -Fi^^2^tan^2/3, 
a2(r7) = — - — sin(2?7)-!^ ^ ^ ^ . (283) 
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D Radiation solutions {k = 0) 

The behaviour of the k = radiation models for this choice of the couphng function is 
determined by 



In 



1 - VU2 



A2/?ln 



Kt] 



r] + 2r]o 



(284) 



where we have exploited our freedom in rjQ to set A = 2Tr]o and U2 is as defined by Eq. (252). 
We also require K=l, so that (f> ^ 4)0 as rj ^ oo. Using Eq. ( ^ ) we can calculate the exact 
evolution of these models 



1 + 



1 + 



V 



-2/13 



2/(3 



(285) 
(286) 



1 Late-time behaviour 

At late times these equations may be approximated by 



a{r]) 



1 - Ai/34 



-7] 



1 



(287) 
(288) 



7] 2 V " / r/^ 

It is necessary to extend the computation of 0(77) to second-order since the first-order 



contributions will later vanish. Eq. ( 287 ) allows asymptotic calculation of conformal time 
as a function of cosmic time, 



(289) 



and hence there is power-law approach to the GR solution 



cP{t) ^ (^0 



2y^2 



/ r \ 1/4 



8 V (in V ) t 



(290) 
(291) 



as t — > 00. 



2 Early-time behaviour 

We now probe the early-time behaviour of these models, finding two distinct cases. 
There exist zeros of as 77 — > and t] — 2t/o. Since we require (p > during the portion 
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of the evolution in which we are interested, and ultimately at late times, we shall take our 
early-time limit to be the most recent of these zeros. We know from Eq. ([TtI) that 



4>o 



1/2 



(292) 



and since (j) € (0,i;^o) and (j) —>■ (j)o at late times we deduce that c/)^ > 0. Prom Eq. ( 292 ), 
we derive X2A > 0. We have fixed A = 2Tr]Q, with F > 0, so we may also say A2??o > 0, 
ie when A2 > 0, r/ ^ is the more recent zero of (p and when A2 < 0, ^ —^Vo is most 
recent. When A2 > we find 



A 



A2 



(2%)"^ V'-'' 



(293) 
(294) 



as 7/ — > 0. We can model the t-parametrised early-time behaviour by Eqs. ( |271| ) - ( |274| ), for 
the k = —1 vacuum models, via the transformation 



(2%) 



A2 



When A2 < 



4^/^<^o 
A 



4V/3</,o 



V + 2r?o 
-2r?o 

' -2r?o)-'^ (r? + 27?o)'+"^ 



(295) 

(296) 
(297) 



as r/ — > —2r]Q. Note that both A and r/o are negative here. The behaviour as a function of t 
can be gleaned from Eqs. (271) - (274), by the redefinitions 



A-2 



7] 

A 



Differentiating Eq. ( ^851) yields 



(f^o (-2r/o) 

-A. 
3 Minima 



2% (1 + A2) + 2r] 



(298) 
(299) 
(300) 



4A27?o 



1 + {v+2vo) 



A2/3 



V r]{r] + 2r]o) 

If the early-time limit is given by X ^ where X is either rj or r] + 2r]Q , then 



(301) 



(302) 
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and thus never vanishes at early times. From Eq. ( p92|) we deduce that when < (p < (pQ, 
(j) is bound to evolve monotonically. The necessity that it tend to a constant at late times 
then occludes the possibility of it diverging after the early-time regime. The non-zero right- 
hand side of Eq. ( p92D on the interval Q < (j) < then guarantees that a is non-zero. 
Consequently, the stationary points in the evolution are obtained by setting the square- 
bracketed factor in Eq. (301) to zero. They lie at rj^, where 

2A2??o 



(r,,+2r;o) 



A2/3 



r]* = r/o (1 + A2) , 



+ 1 



which must be solved numerically for the particular values of /?, A2, and ??o- 



(303) 



E Radiation solutions {k = —1) 

When the curvature is negative, the evolution is determined by the equation 



1 - VU2 



rtanh(7/ -hr/o) + {A^ -T 



A 



rtanh(r/ + 7?o) + {A^ - ^^f'"^ + A 



A2/3 



(304) 



with U2 and A2 as defined earlier. We make the simplifying choice of 

r 

tanh(2r?o) = ^ , 

which ensures that the early-time behaviour (ie, (j) ^ occurs at 7/ = 0. To recover GR at 
late times, we fix 



(305) 



K 



r+(^2_p2)l/2^^ 



r + (^2_r2)i/2_^ • 

Substituting Eq. ( pOS]) into the k = —1 right-hand side of Eq. ( ^) we find 

r A 

yiv) = [cosh(27?) - 1] -I- — sinh(2ry) . 



(306) 



(307) 

Lastly, we note that with r/o given by Eq. ( pOSj ), the expression within the moduli on 
the right-hand side of Eq. ( |304| ) is monotonically increasing from zero and thus positive. 
Dropping the moduli, we find 

^ rtanh(7?+7?o)+(A2-r2)^/^-A \ 
rtanh(??+??o)+(A2-r2)l''2+A I 



1 + K^2l3 ( rtanh(rj+r;o)+(^2_r2)i/2^A 



rtanh(r)+r)o)+(A2-r2)^''-'+yl 



(308) 



4-l//3i^-A2 

2fo 



[r [cosh(2?7) - 1] +Asmh{2r])] 

\ A2/3' 



1 , ^A2/3 ( rtanh(>,+^o)+(A2-r2)^/^-A Y 
rtanh(,?+r?o)+{A2-r2)l''2+A J 



2//3 



rtanh(r;+T;o)+(A2-r2)l/2_^ 
rtanh(r;+r;o)+(A2~r2)i/2+A 



A2 



(309) 
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1 Late-time behaviour 



As ?7 ^ oo we find 



^ XlA^(3 / A-T-{A^-T 



/2 



-4r) 



T + A 

Ho 



,2r) 



1 



2T 



r + A 



-2ri 



(310) 
(311) 



as ^ oo. These lead to the t-dependent behaviours: 

r 



a{t) 



In 



T + A 



Ho 



-t~ 



■^2^V2^ 



_ A|A2^ M + r \ M-r-(A2_r2)^ 
1602 U-rj l^^ + r-(A2_r2)i/2j 

r 



(312) 

(313) 
(314) 



as t — > oo. Again, we see the approacli of this model to the Milne universe as t increases. 

2 Early-time behaviour 

Analysing the behaviour in the neighbourhood of the last zero of (ie, around = 0) 
we obtain a description of the early-time behaviour. We find 



(/»(r/) 



/\r^/PAK~^^ 



A2 



r M2 A M2 



(315) 
(316) 



with K as given by Eq. ( 306] ). The reality condition > and the inequality > 
(from Eq. ( |292| ), using c/;^ > 0) ensure all of the pre- factors in the above relations are 
real and positive. We obtain the t-parametrised behaviour at early times by applying the 
substitutions 



I A, 



A2, 



r M2 



(317) 
(318) 



to Eqs. (|7lD - (^7§ 
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3 Minima 



If A2 > 1 then the scale-factor diverges at early times. For such models to look like GR 
as t — > 00 requires the presence of a minimum. Using Eq. (p7|) we find 



cosh [2(ry^, + ?7o)] = A2 cosh(2r/o) < 



1 - 



1 + K^2/3 



rtanh(T;.+77o)+(yl2- 




-A 


rtanh(r),+»?o)+(-42- 


_r2)i/2_ 




rta.nh{'ij,+rio)+{A'''~ 


_r2)i/2_ 


-A 


_rtanh(r;.+)?o)+(A2- 


_r2)i/2_ 





A2/3 



A2/3 



(319) 



defines r/,,, where the minimum is situated. Although this expression is not soluble analyti- 
cally, we may obtain bounds on the value of 77*. Eq. (p7|), for Theory 3, may be written 



rsinh(277^<) +^cosh(2r/*) = X2AU2 
The bounds on U2, namely U2 G (0, 1) then imply 

< cosh [2 + ??o)] < A2 cosh (2r/o 



1/2 



(320) 



(321) 

after using Eq. ( |305| ). The cosh function is greater than unity and monotonically increasing 
when its argument is positive, allowing us to strengthen the above inequality to 



< V* < -arccosh [A2 cosh (2??o)] — Vo , 
where the lower limit follows from the positivity Of ij. 

F Radiation solutions (k = +1) 

For the positively-curved models {k = +1) we obtain 

F tan(?7 + %) + (r^ + ^2)1/2 -A^^'^ 



(322) 



1 - \/U2 



K 



1 + ^ " Ftan(?7 + ??o) + (r2 + ^2)1/2 + ^ 
This leads to the conformal time-parametrised set of equations 



a{ri) 



(323) 



(t){ri) 



4cj(r/) 



1 




1//3 . 


2410 


Mv) 





r+(r^ + A^y^\m[2{7j + rjo)] 



(324) 
(325) 



G Perfect-fluid solutions (k = 0) 



We examine the late-time form of the solutions resulting from a scalar-tensor cosmology 
driven by the third type of coupling function using the theory first studied by Barrow and 
Mimoso and defined by 

m = ^oew{P^'^) , (326) 
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where P and Q are constants. This choice of arises from the generating function 

1 



and results in the scale-factor 



^3(2-7) 



3(2-7) 
^0 



exp ( -P^' 



(327) 



(328) 



Positivity of the left-hand side of Eq. (p2q) requires that i^PQ > 0. Hence, we may derive 

1 

PQ^ 

and 

4-37 



(329) 



2a; (0 + 3 



PQ ? "I" 2 ? 



3(2-7)2 



J_t2-Q I 4-37 ^2 



(330) 



Again, at late times we demand (p ^ This can be realised as ^ ^ when Q > or as 
^ — > oo when Q < 0. When P = and/or Q = we have GR. 



1 Late-time behaviour 

The limiting forms of the coupling function as — > when Q > 0, and as — > cxd when 
Q < 0, are identical: 



^'^^^^ + '^ 3(2-7)2 pg ^ ' 



In this case, Eq. (32) is 



dt 



2-Q 
2-7 



3(7-1) 
^0 



3(PQ0o)^ 



4_3^\i/2 47-4-Q7 

e exp 



7-1 
2-7 



(331) 



de, (332) 



as ^ 0. This expression may be integrated approximately for small giving 

'7-r 



, , 1 7(2-Q) 

i(0 - 5-^^2(2-^) 



where S is defined by 



1 -P 



7-1 



2-7 



S = 



3j{PQcl)o)— f PQ 



3(7-1) V4-37 



1/2 



(333) 



(334) 



Ignoring the square brackets in Eq. (333), we obtain the lowest-order inversion of this 
expression 

2(2-7) 2(2-7) 

C(t) ~ S'^(2-'3)t^(2-<3) . (335) 

Introducing the corrections associated with the terms in the square bracket leads, at the 
next-order, to 



2(2-7) 2(2-7) 
£,{t) ~ S'^(2-«) fT(2-Q) 



1 + 



2P / 7 - 1 



7 V2-Q 



2Q(2-7) 2Q(2-7) 
5" 7(2-Q) t 7(2-Q) 



(336) 
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From this we can see that ^ — > as t ^ oo iff Q > 2. As ^ — > oo, t ^ oo iff Q < 2; the 
condition </> — > 0o as ^ — > oo requires Q < 0, we thus deduce that the range of theories 
demarcated by < Q < 2 wih not approach GR as t — > cxd. Substituting the above ^{t) 
into Eqs. (p26D and (|32^) we obtain the t-dependent evolution 



<P{t) 



2Q(2-7) 2Q(2-7) 
1 _ PS 7(2-Q) t 7(2-Q) 



ao 



52/37 



.^2/37 



I 2Q(2-7) 2Q(2-7) 

1 PS 7(2-Q) t 7(2-Q) 

3 



(337) 



(338) 



as t 

by 



{PQct>o)'^ 

00. The corresponding late-time evolution of the coupling as a function of (p is given 



2uj{(p) + 3 



4-37 i2-Qf 



(339) 



3(2-7)2 Q In (0/00) 

When Q > the requirement that a be positive ensures P^Q > 0. From Eq. (P) we see 
that for Q > 2 d£,/dt < 0, ie ^ ^ from above and hence P > 0. Eq. ( ^261) confirms 
that (f) ^ (f)Q from above when P and Q lie in these domains. Conversely, when Q < the 
requirement a > implies P^'^ < 0. At late times ^ — > 00 in these models and so P < 
and 4> approaches the GR value, 0o, from below. 

We remarked earlier that the Q = case gives pure GR at all times. When Q = 2 
Eq. ( 330 ) for 2a; + 3 does not approach the limit given in Eq. ( |331|) as ^ — > 0, instead we 
have 

2(4 - 37)3p 



2w(0 + 3 



(340) 



12(2-7)2 

and 2a; + 3 decays to zero as ^ — > 0. The range of Q incompatible with current observations, 
that GR is a good approximation to the time description of gravitation today, is then 
< Q < 2. 



2 Early-time behaviour 

When Q > 2, P > the scale-factor a approaches zero as — > ±00. 



ie, BD theory. The square root of this limit is positive and so it follows from Eq. ( |32| ) that 
d^/dt must also be positive, demanding that ^ approach —00. The positivity of a, Q and 
P imply (—1)*^ > and hence ^ 00 as ^ — > — cxd, a ^ 0. 

When Q < 0, a necessary condition for the scale-factor to approach zero is ^ ^ 
and in this limit we recover BD theory again, the coupling given by Eq. ( [3411) . As before 



\/2a7+~3 > 0, d^^/dt > and ^ from above. Positivity of a thus requires PQ > 0, 
which is guaranteed from the late-time behaviour. However, for a to converge to zero as 
,^ ^ requires P > and hence Q > 0. This direct contradiction ensures that a(^) will be 
non-singular when Q < 0. Explicitly, a{^) possesses a minimum at 
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For reasons given in Section we refrain from presenting the explicit form of the solution 
here. 



3 Dust models 

The late-time evolution of the 7 = 1 models is, from Eqs. ( p37D and (338), that 



ait) 



1 - p 



( 3 

V2(PQ)V2 



2Q 



«0 



/ 9 



1/3 



^2/3 



p 

1 — 



3 V2(PQ)i/2 



2Q 

2^ 2Q 
t2-Q 



(343) 
(344) 



as t — > 00. 



4 Inflationary models 

As was the case in Section , the solutions when 7 = are qualitatively different to their 
more general counterparts, Eqs. ( p37| ) and ( |33§| ). In this case Eq. ( p32[ ) is 



dt ~ 



1/2 

exp 



2^ 



as — > 0, which integrates approximately in this limit to give 



(2 



3a|] 



1/2 

^^ine 



1 + 



2gin^ 



The first-order inversion of this at large t gives 

^(t) ~ exp 



3alt 



,(2 



1/2 



3a[]t 



.(2 



1/2 



3Qa[]t 



1/2 
- 



as t ^ cxo. Substituting this into Eqs. (|326[) and ( |328|) yields 

SQa^t 



ait) 



1 + P exp 



.(2 



1/2 



(345) 



(346) 



(347) 



which arises by neglecting the square brackets on the right of Eq. ( |346| ). The next-order 
correction to this result is 



(348) 



ao 



a^t 



, Pi2 + Q) 
' ^^"P 



3Qalt 



1/2 



(349) 
(350) 



51 



5 Connection to the parameters of Theory 3 

At late times we have for all of the models considered in this section 



2u!{(p) + 3 



4-37 {2-Q)^E 
3(2-7)2 Q 

4-37 i2-Q)^E 
3(2-7)2 i-Q) 



1 -1 



1 



, as 



(351) 
(352) 



where E is a constant. When (p > (po and Q > and we can model the form of the coupling 
as a function for Theory 3 by Eq. ( |35lD . When cj) < 4>o and Q < we can use Eq. (|352D 
for the coupling as a function of (p. This leads to the consistency relations 



(2 - Qf 
sgn(Q) 



E 



3(2-7)2 53 
4-37 /? ' 



connecting Theory 3 as defined in Section IV with the solutions presented here. 



6 Exact solution 

We note the existence of an exact solution, defined by the choice 

m = <Ao (i - e-^) , 

where A is a constant. This leads, by Eq. (|4^), to 

1 



9(0 



A-2 



and by Eq. (|45|) to 



3(2-7^ 



,3(2-7) ^ "O ) t\ 



(A - 2)^0 

The exact form of coupling driving this behaviour is given by Eq. ( ^0|) as 

2 



2c^(e) +3 



4-37 



A^^-2 _ 2j 



3(2-7)2(A-2) 



or, using Eq. (p5l 
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4-37 
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where 



3(2-7)2(A-2) (i_^/0o)-i_ J 
J = l I (37-4)(A-2) 



(353) 
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VIII Discussion 



In this paper we have supphed a comprehensive study of isotropic cosmological models in 
scalar-tensor theories, extending the earlier work of refs.Jl^ and |27|. We have explored the 
behaviour of isotropic cosmological models in these theories using a combination of two basic 
mathematical techniques introduced in section ED . In the case where the trace of the energy- 
momentum tensor of matter vanishes (ie vacuum or radiation) exact solutions can be found 
directly for all curvatures if the requisite integrals can be performed, asymptotic forms are 
easily derived in all cases. This technique exploits the conformal relationship between scalar- 
tensor theories and general relativity that exists when the trace of the energy-momentum 
tensor vanishes. However, when the energy-momentum tensor is not trace-free the conformal 
equivalence disappears and the indirect method of Barrow and Mimoso must be used to 
find exact solutions. This only works for zero curvature cosmological models but includes 
the important cases of all p = universes and inflationary universes with —p < p < — 
It also permits a simple means of comparing the behaviour of cosmological solutions in any 
scalar-tensor theory with those of Brans-Dicke theory at early and late times. Since this 
procedure does not commence with the specification of uj{(j)), but with a choice of generating 
function that produces the entire solution by non-linear transformations, it is necessary to 
build up intuition by a thorough exploration of the results of employing particular classes 
of generating function. In particular, we were able to find generating functions which gave 
rise to dust universes in scalar-tensor theories for which the exact radiation and vacuum 
solutions can be found exactly by the direct method. This provides us with descriptions of 
scalar-tensor cosmologies throughout the entire radiation and dust and vacuum dominated 
eras. We were also able to find a wide range of new inflationary universe solutions with 
p = —p in these theories. 

In section [V we introduced three classes of scalar-tensor theory which permit asymptotic 
approach to general relativity at late times when cj) ^ <j)Q. The parameters defining the 
functional form of w((/>) which specify these gravity theories can be restricted further if 
we require the theory to approach general relativity in the weak-field limit (lj ^ oo and 
u}'uj~^ — > 0), and describe expanding universes. For each of these general classes of theory 
we have determined the behaviour of flat, open, and closed universes by a combination of 
exact solutions and asymptotic studies of the early and late-time behaviours. 

The catalogue of solutions and asymptotes that we have found will enable scalar-tensor 
theories to be constrained in new ways because they enable complete cosmological histo- 
ries to be constructed through initial vacuum, radiation, dust and final vacuum-dominated 
eras. The standard sequence of physical processes responsible for events like monopole pro- 
duction, inflation, baryosynthesis, primordial black hole formation, electroweak unification, 
the quark-hadron phase transition, and nucleosynthesis can be explored in the cosmologi- 
cal environment provided by scalar-tensor gravity theories. The constraints derived from 
these considerations can be compared directly with those imposed by weak-field tests in 
the solar system and observations of astrophysical objects like white dwarfs and the binary 
pulsar. The ubiquity of scalar fields in current string theories of high-energy physics has led 
to continued interest in the detailed behaviour of scalar-tensor gravity theories and their 
associated cosmologies. In this paper we have displayed some of the diversity that these 
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cosmologies possess together with a collection of methods for solving other specific theories 
that may be motivated by future developments in high-energy physics. 
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Appendix — Conformal equivalence 

The action for trace-free matter is conformally invariant and we may exploit this fact 
to study the behaviour of vacuum and radiation models with (p > (f)Q. Under a conformal 
transformation to a new metric g, with components 



5ab 



(,-2 



and a re-definition of the field 



the action becomes 



Sg 



d x\J —g 



(361) 
(362) 

(363) 



neglecting overall constant factors. When < < i;^o the coupling function for Theory 2, 
as defined in section IV, is 



2w((^) + 3 = ^2 



and so 



2cj( 



+ 3 = ^2 



ln( ^ 

30 



In 



-25 



-25 



(364) 



(365) 



with (j) > (pQ when < (j) < (pQ. Eq. (365) is exactly the form of the coupling for Theory 
2 with a gravitational scalar field (j) > (j)Q. Thus we may obtain solutions when (j) > (p^ m. 
Theory 2 by applying the transformations 

o — > 4>a , (366) 
<t> ^ \, (367) 



in this order, to the solutions presented in Section |V|. These are transformations which 



render the form of y in Eq. (15) invariant. 

Lastly, we remark that the asymptotic behaviours of Theories 1 and 3 when cp > (pQ may 
also be examined in this way, since the coupling functions for both these theories may be 
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approximated by logarithms as — > 0o- 
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